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Abstract 

In an atomic gas near a Feshbach resonance, the energy of two colliding atoms is 
close to the energy of a bound state, i.e., a molecular state, in a closed channel that 
is coupled to the incoming open channel. Due to the different spin arrangements 
of the atoms in the open channel and the atoms in the molecular state, the energy 
difference between the bound state and the two-atom continuum threshold is exper- 
imentally accessible by means of the Zeeman interaction of the atomic spins with a 
magnetic field. As a result, it is in principle possible to vary the scattering length to 
any value by tuning the magnetic field. This level of experimental control has opened 
the road for many beautiful experiments, which recently led to the demonstration 
of coherence between atoms and molecules. This is achieved by observing coherent 
oscillations between atoms and molecules, analogous to coherent Rabi oscillations 
that occur in ordinary two-level systems. We review the many-body theory that 
describes coherence between atoms and molecules in terms of an effective quan- 
tum field theory for Feshbach-resonant interactions. The most important feature of 
this effective quantum field theory is that it incorporates the two-atom physics of 
the Feshbach resonance exactly, which turns out to be necessary to fully explain 
experiments with Bose-Einstein condensed atomic gases. 
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1 Introduction 



Following the first experimental realization of Bose-Einstein condensation [1], 
a great deal of experimental and theoretical progress has been made in the 
field of ultracold atomic gases [2,3,4,5]. One particular reason for this progress 
is the unprecedented experimental control over the atomic gases of interest. 
This experimental control over the ultracold magnetically-trapped alkali gases, 
has recently culminated in the demonstration of experimentally adjustable 
interactions between the atoms [6]. This is achieved by means of a so-called 
Feshbach resonance [7]. 

Feshbach resonances were introduced in nuclear physics to describe the nar- 
row resonances observed in the total cross section for a neutron scattering of 
a nucleus [8]. These very narrow resonances are the result of the formation of 
a long-lived compound nucleus during the scattering process, with a binding 
energy close to that of the incoming neutron. The defining feature of a Fesh- 
bach resonance is that the bound state responsible for the resonance exists in 
another part of the quantum-mechanical Hilbert space than the part associ- 
ated with the incoming particles. In the simplest case, these two parts of the 
Hilbert space are referred to as the closed and open channel, respectively. 

Following these ideas from nuclear physics, Stwalley [9] and Tiesinga et al. [10] 
considered Feshbach resonances in ultracold doubly spin-polarized alkali gases. 
Due to the low temperatures of these gases, their effective interatomic inter- 
actions are to a large extent completely determined by the s-wave scattering 
length. Analogous to the formation of a compound nucleus in neutron scat- 
tering, two atoms can form a long-lived bound state, i.e., a diatomic molecule, 
during an s-wave collision. This process is illustrated in Fig. 1. The two incom- 
ing atoms in the open channel have a different hyperfine state than the bound 
state in the closed channel and the coupling between the open and closed 
channel is provided by the exchange interaction. As a result of this difference 
in the hyperfine state, the two channels have a different Zeeman shift in a 
magnetic field. Therefore, the energy difference between the closed-channel 
bound state and the two-atom continuum threshold, the so-called detuning, is 
experimentally adjustable by tuning the magnetic field. This implies that the 
s-wave scattering length, and hence the magnitude and sign of the interatomic 
interactions, is also adjustable to any desirable value. In Fig. 2 the scattering 
length, as measured by Inouye et al. [6], is shown as a function of the mag- 
netic field. The position of the resonance in the magnetic field is at B ~ 907 
(G)auss in this case. Following this first experimental observation of Feshbach 
resonances in 23 Na [6], they have now been observed in various bosonic atomic 
species [11,12,13,14,15], as well as a number of fermionic isotopes [16,17,18,19]. 

With this experimental degree of freedom it is possible to study very interest- 
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Fig. 1. Illustration of a Feshbach-resonant atomic collision. Two atoms, with a 
hyperfine state indicated by the arrow, collide and form a long-lived molecule with 
a different spin arrangement, which ultimately decays again into two atoms. 

ing new regimes in the many-body physics of ultracold atomic gases. The first 
experimental application was the detailed study of the collapse of a condensate 
with attractive interactions, corresponding to negative scattering lengths. In 
general a collapse occurs when the attractive interactions overcome the sta- 
bilizing kinetic energy of the condensate atoms in the trap. Since the typical 
interaction energy is proportional to the density, there is a certain maximum 
number of atoms above which the condensate is unstable [20,21,22,23,24]. In 
the first observations of the condensate collapse by Bradley et al. [25] , a con- 
densate of doubly spin-polarized 7 Li atoms was used. In these experiments the 
atoms have a fixed negative scattering length which for the experimental trap 
parameters lead to a maximum number of condensate atoms that was so small 
that nondestructive imaging of the condensate was impossible. Moreover, ther- 
mal fluctuations due to a large thermal component made the initiation of the 
collapse a stochastic process [26], thus preventing also a series of destructive 
measurements of a single collapse event [27]. A statistical analysis has never- 
theless resulted in important information about the collapse process [28] . Very 
recently, it was even possible to overcome these complications [29]. 

In addition to the experiment with 7 Li, experiments with 85 Rb have been 
carried out [30]. In particular, Roberts et al. [31] also studied the stability 
criterion for the condensate, and Donley et al. [32] studied the dynamics of 
a single collapse event in great detail. Both of these experiments make use 
of a Feshbach resonance to achieve a well-defined initial condition for each 
destructive measurement. It turns out that during a collapse a significant 
fraction of atoms is expelled from the condensate. Moreover, one observes 
a burst of hot atoms with an energy of about 150 nK. Several mean-field 
analyses of the collapse, which model the atom loss phenomenologically by a 
three-body recombination rate constant [33,34,35,36,37,38,39], as well as an 
approach that considers elastic condensate collisions [40,41], and an approach 
that takes into account the formation of molecules [42], have offered a great 
deal of theoretical insight. Nevertheless, the physical mechanism responsible 
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Fig. 2. The scattering length as a function of magnetic field as measured by Inouye 
et al. [6]. The scattering length is normalized such that it is equal to one far off 
resonance. 

for the explosion of atoms out of the condensate and the formation of the 
noncondensed component is to a great extent still not understood at present. 

A second experimental application of a Feshbach resonance in a Bose-Einstein 
condensed gas is the observation of a bright soliton train by Strecker et al. 
[15]. In this experiment, one starts with a large one-dimensional Bose-Einstein 
condensate of 7 Li atoms with positive scattering length near a Feshbach res- 
onance. The scattering length is then abruptly changed to a negative value. 
Due to its one-dimensional nature the condensate does not collapse, but in- 
stead forms a train of on average four bright solitary waves that repel each 
other. The formation of these bright solitons is the result of phase fluctu- 
ations [43], which are in this case important due to the low dimensionality 
[44,45,46,47,48,49]. The repulsion between the bright solitons is a result of 
their relative phase difference of about ir. In a similar experiment Khaykovich 
et al. [50] have observed the formation of a single bright soliton. 

A third experimental application are the experiments with trapped gases of 
fermionic atoms, where the objective is to cool the gas down to temperatures 
where the so-called BCS transition, i.e., the Bose-Einstein condensation of 
Cooper pairs, may be observed. The BCS transition temperature increases if 
the scattering length is more negative [51], and hence a Feshbach resonance can 
possibly be used to make the transition experimentally less difficult to achieve. 
This possibility has inspired the study of many-body effects in fermionic gases 
near a Feshbach resonance [52,53,54,55,56,57,58], as well as fluctuation effects 
on the critical temperature [59,60]. One of the most interesting features of a 
fermionic gas near a Feshbach resonance is the crossover between a conden- 
sate of Cooper pairs and a condensate of molecules, the so-called BCS-BEC 
crossover that was recently studied by Ohashi and Griffin [55,56,57] on the 
basis of the Nozieres-Schmitt-Rink formalism [61]. As a first step towards this 
crossover, Regal et al. [62] were recently able to convert a fraction of the atoms 
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in a gas of fermionic atoms in the normal state into diatomic molecules, by 
sweeping the magnetic field across a Feshbach resonance. Following this obser- 
vation, Strecker et al. observed the formation of long-lived 6 Li 2 molecules [63], 
and Xu et al. observed 23 Na molecules [64]. Very recently, even the formation 
of Bose-Einstein condensates of molecules has been observed by Jochim et al. 
[65], Greiner et al. [66], and by Zwierlein et al. [67]. As another application of 
Feshbach resonances in fermionic gases we mention here also the theoretical 
proposal by Falco et al. to observe a new manifestation of the Kondo effect in 
these systems [68]. 

The experimental application on which we focus in this paper is the obser- 
vation of coherent atom-molecule oscillations [69]. These experiments are in- 
spired by the theoretical proposal of Drummond et al. [70] and Timmermans et 
al. [71] to describe the Feshbach-resonant part of the interactions between the 
atoms in a Bose-Einstein condensate by a coupling of the atomic condensate to 
a molecular condensate. For this physical picture to be valid, there has to be a 
well-defined phase between the wave function that describes the atoms in the 
atomic condensate, and its molecular counterpart. An equivalent statement is 
that there is coherence between the atoms and the molecules. Since the energy 
difference between the atoms and the molecular state is experimentally tun- 
able by adjusting the magnetic field, it is, with this physical picture in mind, 
natural to perform a Rabi experiment by means of one pulse in the magnetic 
field towards resonance, and to perform a Ramsey experiment consisting of 
two short pulses in the magnetic field. If the physical picture is correct we 
expect to observe oscillations in the remaining number of condensate atoms 
in both cases. 

In the first experiment along these lines, Claussen et al. [72] started from a 
Bose-Einstein condensate of 85 Rb atoms without a visible thermal cloud and 
tuned the magnetic field such that the atoms were effectively noninteracting. 
With this atomic species this is possible, because the off-resonant background 
scattering length is negative, which can be compensated for by making the res- 
onant part of the scattering length positive. Next, one applied a trapezoidal 
pulse in the magnetic field, directed towards resonance. As a function of the 
duration of the pulse one observed that the number of atoms first decreases 
but after some time increases again. This increase can not be explained by a 
"conventional" loss process, such as dipolar relaxation or three-body recombi- 
nation, since the magnitude of the loss is in these cases given by a rate constant 
times the square and the cube of the density, respectively. As a result, the loss 
always increases with longer times. A theoretical description of this experi- 
ment is complicated by the fact that the experiment is at long times close to 
the resonance where little is known about the magnetic-field dependence of 
these rate constants. Although the magnetic-field dependence has been calcu- 
lated for a shape resonance [73,74,75,76], it is not immediately obvious that 
the results carry over to the multichannel situation of a Feshbach resonance. 
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Moreover, precise experimental data is unavailable [77]. Therefore a satisfy- 
ing quantitative description is still lacking, although two attempts have been 
made [41,78]. 

After these experiments, the same group performed an experiment consisting 
of two short pulses in the magnetic field towards resonance, separated by a 
longer evolution time [69]. As a function of this evolution time an oscillation 
in the number of condensate atoms was observed. Over the investigated range 
of magnetic field during the evolution time, the frequency of this oscillation 
agreed exactly with the molecular binding energy found from a two-atom 
coupled-channels calculation [79], indicating coherence between atoms and 
molecules. Very recently, Claussen et al. have performed a similar series of 
measurements over a larger range of magnetic fields [80]. It was found that 
close to resonance the frequency of the oscillation deviates from the vacuum 
molecular binding energy as a result of many-body effects [81,82]. 

As already mentioned, the first theories for Feshbach-resonant interactions in- 
troduce the physical picture of an interacting atomic Bose- Einstein condensate 
coupled to a noninteracting molecular condensate [70,71,83]. The first descrip- 
tion of the Ramsey experiments by Donley et al. [69] was achieved within the 
Hartree-Fock-Bogoliubov mean-field theory [79,78,84]. 

It turns out that, for a complete understanding of the experiments, it is neces- 
sary to exactly incorporate the two-atom physics into the theory. Although the 
above-mentioned theories have provided a first understanding of the physics of 
a Bose gas near a Feshbach resonance, these many-body theories do not con- 
tain the two-atom collision properties exactly. To incorporate the two-atom 
physics exactly, it is from a diagrammatic point of view required to sum all 
the ladder Feynman diagrams of the microscopic theory. By means of this 
procedure, we have recently derived an effective quantum field theory describ- 
ing the many-body properties of an atomic gas near a Feshbach resonance 
[85]. It is the aim of this paper to review and extend this effective atom- 
molecule theory and its applications [85,81,82]. Moreover, along the way we 
discuss some of the differences and similarities between our theory and a num- 
ber of other theories for Feshbach-resonant interactions in atomic Bose gases 
[70,71,83,79,78,84,86,87,88,89,90]. 

With this objective in mind, this paper is organized as follows. In Section 2 we 
review two-atom scattering theory. In particular, we emphasize the relation 
between the scattering amplitude of a potential and its bound states. Both 
the single-channel well as the multichannel case that can give rise 

to Feshbach resonances, are discussed. This introductory section introduces 
many important concepts in a simple setting, and hence clarifies much of the 
physics that is discussed in later sections. In Section 3 we present in detail 
the derivation of an effective quantum field theory applicable for studying 
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many-body properties of the system, starting from the microscopic atomic 
hamiltonian for a Feshbach resonance. This effective field theory consists of an 
atomic quantum field that is coupled to a molecular quantum field responsible 
for the Feshbach resonance. It is used in Section 4 to study the normal state of 
the gas. In particular, we show here that the two-atom scattering properties as 
well as the molecular binding energy are correctly incorporated into the theory. 
Moreover, we also discuss many-body effects on the molecular binding energy. 
Section 5 is devoted to the discussion of the Bose-Einstein condensed phase 
of the gas. We derive the mean-field theory resulting from our quantum field 
theory. We also discuss the differences and similarities between this mean-field 
theory and in particular the mean-field theories that were recently proposed 
by Kokkelmans and Holland [79], Mackie et al. [78], and Kohler et al. [84]. 
In Section 6 our mean-field theory is applied to the two-pulse experiments 
[69,80]. It is the perfect agreement between theory and experiment obtained 
in this section that ultimately justifies the ah initio approach to Bose gases 
near a Feshbach resonance reviewed in this paper. We end in Section 7 with 
our conclusions. 



2 Scattering and bound states 

In this section we give a review of quantum-mechanical scattering theory. We 
focus on the relation between the scattering amplitude of a potential and its 
bound states [91,92]. In the first part we consider single-channel scattering and 
focus on the example of the square well. In the second part we consider the 
situation of two coupled channels, which can give rise to a Feshbach resonance. 

2.1 Single- channel scattering: an example 

We consider the situation of two atoms of mass m that interact via the po- 
tential V(r) that vanishes for large distances between the atoms. The motion 
of the atoms separates into the trivial center-of-mass motion and the relative 
motion, described by the wave function ip(r) where r = x x — x 2 , and x x and 
x 2 are the coordinates of the two atoms, respectively. This wave function is 
determined by the time-independent Schrodinger equation 

V>(r) = E^(r) , (1) 

with E the energy of the atoms in the center-of-mass system. Solutions of 
the Schrodinger equation with negative energy correspond to bound states of 
the potential, i.e., to molecular states. To describe atom-atom scattering we 
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Fig. 3. Schematic representation of two-atom scattering in the center-of-mass refer- 
ence frame. The atoms are initially in a plane-wave state with relative momentum 
Kk, and scatter into the spherical wave with relative momentum hk' . Due to energy 
conservation we have that k = k' . The angle between k and k' is denoted by 6. The 
region where the interaction takes place is indicated by the black circle. 



have to look for solutions with positive energy E = 2ek, with = h 2 k 2 /2m 
the kinetic energy of a single atom with momentum hk. Since any realistic 
interatomic interaction potential vanishes rapidly as the distance between the 
atoms becomes large, we know that the solution for r — > oo of Eq. (1) is given 
by a superposition of incoming and outgoing plane waves. More precisely, the 
scattering wave function is given by an incoming plane wave and an outgoing 
spherical wave and reads 



Ak'r 



^(r)~e ik ' r + /(k',k)- 



(2) 



where the function /(k',k) is known as the scattering amplitude. The inter- 
atomic interaction potential depends only on the distance between the atoms 
and hence the scattering amplitude depends only on the angle 6 between k 
and k' = k'r, and the magnitude k. Because of energy conservation we have 
that k! = k. The situation is shown schematically in Fig. 3. 

Following the partial-wave method we expand the scattering amplitude in 
Legendre polynomials Pi(x) according to 



/(tf,k) = £/,(*)fl(cos0) . 

1=0 



(3) 



The wave function is expanded in a similar manner as 



il>(r,e) = '£R l (k,r)P l (cose) , 
1=0 



(4) 



with Ri(k,r) = ui(k,r)/r the radial wave function and ui(k,r) determined by 
the radial Schrodinger equation 



d 2 1(1 + 1) mV{r) 



dr 2 



h 2 



+ k 2 



ui(k,r) = . 



(5) 



By expanding also the incident plane wave in partial waves according to 



1=0 



kr 



(6) 



we can show that to obey the boundary condition in Eq. (2), the partial- wave 
amplitudes fi(k) have to be of the form 



fi(k) 



21 + 1 

2ik 



(7) 



where 5i(k) is the so-called phase shift of the l-th partial wave. 

For the ultracold alkali atoms, we are allowed to consider only s-wave (/ = 0) 
scattering, since the colliding atoms have too low energies to penetrate the 
centrifugal barrier in the effective hamiltonian in Eq. (5). Moreover, as we 
see later on, the low-energy effective interactions between the atoms are fully 
determined by the s-wave scattering length, defined by 



a = — lim 



So(k) 



klO k 



(8) 



From Eq. (7) we find that the s-wave scattering amplitude is given by 
1 



fo(k) 



k cot S (k) — ik 



(9) 



As explained above, we take only the s-wave contribution into account, which 
gives for the scattering amplitude at zero-momentum 



/(0,0)~-a. 



(10) 
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To illustrate the physical meaning of the s-wave scattering length, we now 
calculate it explicitly for the simple case that the interaction potential is a 
square well. We thus take the interaction potential of the form 



V(r) 



V ifr< R; 
if r > R, 



with R > 0. With this potential, the general solution of Eq. (5) for I 
given by 



is 



u<(r) = Ae ik<r + Be~ ik< ~ r \ for r < R; 



k<T 



u >(r) = Ce ikr + De 



-ikr 



for r > R, 



[12) 



with k K = \Jk 2 — mVo / K 2 . Since the wave function ip{r) has to obey the 
Schrodinger equation at the origin we have to demand that the function ■u < (r) 
vanishes at this point. This leads to the boundary condition B = —A. By com- 
paring the explicit form of the wave function w > (r) with the s-wave component 
of the general scattering wave function for r — > oo, we find that 



3 2i5 (fc) _ _ 



c 

D 



(13) 



Hence, we determine the phase shift by demanding that the wave functions 
for r < R and r > R join smoothly. This leads to the equations 



^ ^ e ik<R _ e -ik<R^j _ _ e 2i8 (k) e ikR _^ & -ikR 

A (k<e lk<R + A;<e" ifc<i? ) = - e ™ (k) ke ikR _ ke ~ikR ^ 



(14) 



where we have chosen the normalization such that D — 1. Multiplication of 
the above equations with e~ lS °^ and dividing the result leads to 



fctan(A; < i?) = k< tan(5 (A;) + kR) , 



(15) 



from which it follows that 



S (k) = -kR + tan' 1 



— tan(fc < i2) 



(16) 



Note that for a repulsive hard-core potential we have that V — > oo and there- 
fore, with the use of the definition in Eq. (8), that the scattering length a = R. 
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Fig. 4. Scattering length (solid line) and effective range (dashed line) for an at- 
tractive square well in units of the range of the potential, as a function of the 
dimensionless parameter 7 = R^m\Vo\/h 2 . 

This immediately gives a physical picture for a positive s-wave scattering 
length: at low energy and momenta the details of the potential are unimpor- 
tant and we are allowed to model the potential with an effective hard-core 
potential of radius a. For a fully repulsive potential the scattering length is 
always positive. For a potential with attractive parts the scattering length 
can be both negative and positive, corresponding to attractive and repulsive 
effective interactions, respectively. 

This is seen by explicitly calculating the scattering length for our example in 
the case that V < 0. As its definition in Eq. (8) shows, the scattering length 
is determined by the linear dependence of the phase shift on the magnitude 
of the relative momentum hk of the scattering atoms for small momentum. 
Generally, the phase shift can be expanded according to [91,92,93] 



kcot(8 (k)) = -- + \r eS k 2 + 
a I 



(17) 



from which the scattering length is determined by 
tan7\ 



a = R[l 



7 



(18) 



with 7 = Rym\V \/h 2 a dimensionless constant. The parameter r eff is the 
so-called effective range and is, in our example of the square-well potential, 
given by 



1 + 



3tan7-7(3 + 7 2 ) 
37(7 — tan7) 2 



(19) 
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In Fig. 4 the scattering length is shown as a function of 7 by the solid line. 
Clearly, the scattering length can be both negative and positive, and becomes 
equal to zero at values of 7 such that 7 = tan 7. In the same figure, the 
effective range is shown by the dashed line. Note that the effective range 
diverges if the scattering length becomes equal to zero. This is because the 
expansion in Eq. (17) is ill-defined for a = 0. At values of 7 = (n + l/2)n 
with n a positive integer the scattering length diverges and changes sign. This 
behaviour is called a potential or shape resonance and in fact occurs each time 
the potential is just deep enough to support a new bound state. Therefore, 
for large and positive scattering length the square well has a bound state 
with an energy just below the continuum threshold. It turns out that there 
is an important relationship between the energy of this bound state and the 
scattering length. 

To find this relation we have to determine the bound-state energy by solv- 
ing the Schodinger equation for negative energy V < E < 0. This leads to 
solutions 

u<(r) = A (e ik<r - e~ ik<r ) , for r < R; 

^ V J ' (20) 

w>(r) = Be~ Kr , for r > R, 



with k K = \jm{E — V )/h 2 and n = ym\E\/h 2 . Demanding again that these 
solutions join smoothly at r — R, we find the equation for the bound-state 
energy 

y/^EJ = -y/%(E m -Vo)cot (J^(E m - V )) . (21) 



We can show that for values of 7 such that (n — l/2)n < 7 < (n + l/2)n this 
equation has n solutions for V < E m < [92]. 

For small binding energy \E m \ <ti \Vq\ we have from the equation for the 
bound-state energy that 



/ Tfi 

— |£ m | ~ - 7 cot7/i?~ l/a , (22) 



where we made use of the fact that 7 has to be close to the resonant values 
(n + 1/2) 7r in this case. This leads to the desired relation between the energy 
of the molecular state and the scattering length given by 
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This result does not depend on the specific details of the potential and it turns 
out to be quite general. Any potential with a large positive scattering length 
has a bound state just below the continuum threshold with energy given by 
Eq. (23). Moreover, the relation will turn out to hold also in the multichannel 
case of a Feshbach resonance as we will see in Section 2.3. Before discussing 
this situation, we first turn to some concepts of scattering theory which are of 
importance for the remainder of this paper. 



2.2 Single- channel scattering: formal treatment 



Let us give a more formal treatment of the scattering theory described above. 
In a basis-independent formulation the Schrodinger equation we have solved 
reads 



H + V 



E\ 



(24) 



with H — P 2 / m the relative kinetic energy operator for the atoms. To describe 
scattering, we have to look for solutions which asymptotically represent an 
incoming plane wave, and an outgoing spherical wave. In the absence of the 
potential V there is no scattering, and hence we demand that the solution 
of Eq. (24) reduces to a plane wave in the limit of vanishing potential. The 
formal solution that obeys this condition is given by 

= |k> + —^V\^ ] ) , (25) 



where |k) represents the incoming plane wave and we recall that E = 2e^ 
is the kinetic energy of the atoms. This energy is made slightly complex by 
the usual limiting procedure E + = lim^o E + irj. Moreover, we have for the 
scattering amplitude that 

/(k / ,k) = -£^(k'|V|4 +) ). (26) 

To determine the scattering amplitude directly, we introduce the two-body 
T(ransition) matrix by means of 

Wi +) >=T 2B (£+)|k> . (27) 

Multiplying the formal solution in Eq. (25) by V we have that 

f 2B (£ + )|k) = V\k) + V 1 -^f 2B (E + )\k) . (28) 

E+ — H 
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Since this equation holds for an arbitrary plane wave |k) and because these 
plane waves form a complete set of states we have the following operator 
equation for the two-body T-matrix 

f 2B (z) = V + V—L-T 2B (z) . (29) 
z-H 



This equation is called the Lippmann-Schwinger equation and from its solution 
we are able to determine the scattering properties of the potential V. To see 
this we first note that from the definition of the T-matrix in Eq. (27), together 
with Eq. (26), it follows immediately that 

/(k/ ' k) = "4^ (k ' |:f2B(24)|k) - (30) 



Therefore, we indeed see that the two-body T-matrix completely determines 
the scattering amplitude. The Lippmann-Schwinger equation for the two-body 
T-matrix can be solved in perturbation theory in the potential. This results 
in the so-called Born series given by 



f 2B {z) = V + VG (z)V + VG (z)VG (z)V + ■■■ , (31) 
where 

G {z) = — , (32) 
z- H n 



is the noninteracting propagator of the atoms. By using, instead of the true 
interatomic interaction potential, a pseudopotential of the form 

V(x-x / ) = ^^5(x-x / ) , (33) 
m 



the first term in the Born series immediately yields the correct result for the 
scattering amplitude at low energies and momenta, given in Eq. (10). Such a 
pseudopotential should therefore not be used to calculate higher-order terms 
in the Born series, but should be used only in first-order perturbation theory. 

The poles of the T-matrix in the complex-energy plane correspond to bound 
states of the potential. To see this we note that the formal solution of the 
Lippmann-Schwinger equation is given by 

f 2B (z) = V + V—^V . (34) 
z — H 
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After insertion of the complete set of eigenstates \ip a ) of H = H + V we have 
f 2B (z) = V + J2 V ^ a) ^ V , (35) 



where the summation over a is discrete for the bound-state energies e a < 0, 
and represents an integration for positive energies that correspond to scatter- 
ing solutions of the Schrodinger equation, so explicitly we have that 



I°*(z) = V + Y, V l -^Mv + / JLy l^ +) )^ + V , (36) 



(2vr) 3 z - 2e k 

From this equation we clearly see that the two-body T-matrix has poles in the 
complex-energy plane, corresponding to the bound states of the potential. In 
addition, the T-matrix contains a branch cut on the positive real axis due to 
the continuum of scattering states. 

As an example, we note that for s-wave scattering the T-matrix T 2B (2e£) = 
(k'|T 2B (2e£)|k) is independent of the angle between k' and k. From the rela- 
tion between the T-matrix and the scattering amplitude, and the expression 
for the latter in terms of the phase shift, we have for low positive energies 



T 2B (E + ) 



Arrh 2 



m 



Anah 2 



m 



1 



1 _l_ in mE ar c gmE 



(37) 



where we made use of the expansion in Eq. (17). From this result we deduce 
by analytic continuation that 



T 2B (z) ~ 



Anah 2 



m 



lyf- 



mz ar c ffinz 

' h 2 2h 2 



(38) 



Clearly, for large and positive scattering length the T-matrix has a pole at 
negative energy E m = — h 2 /ma 2 , in complete agreement with our previous 
discussions. 

Summarizing, we have found that the scattering length of an attractive poten- 
tial well can have any value and depends strongly on the energy of the weakliest 
bound state in the potential. In principle therefore, if we have experimental 
access to the energy difference of this bound state and the continuum thresh- 
old we are able to experimentally alter the scattering length and thereby the 
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effective interactions of the atoms. In the single-channel case this is basically 
impossible to achieve. In a multichannel system, however, the energy difference 
is experimentally accessible, which makes the low-energy effective interactions 
between the atoms tunable. In the next section we discuss this situation. 



2.3 Example of a Feshbach resonance 



We consider now the situation of atom-atom scattering where the atoms have 
two internal states [94]. These states correspond, roughly speaking, to the 
eigenstates of the spin operator S of the valence electron of the alkali atoms. 
The effective interaction potential between the atoms depends on the state 
of the valence electrons of the colliding atoms. If these form a singlet the 
electrons are in principle allowed to be on top of each other. For a triplet this 
is forbidden. Hence, the singlet potential is generally much deeper than the 
triplet potential. 

Of course, in reality the atom also has a nucleus with spin I which interacts 
with the spin of the electron via the hyperfine interaction 

V M = ^1 • S, (39) 



with cihf the hyperfine constant. The hyperfine interaction couples the singlet 
and triplet states. Moreover, in the presence of a magnetic field the different 
internal states of the atoms have a different Zeeman shift. In an experiment 
with magnetically-trapped gases, the energy difference between these states 
is therefore experimentally accessible. Putting these results together, we can 
write down the Schodinger equation that models the above physics 



■*£- + Vt{t)-E Vm 
h 2 v 2 



V"r(r) 



V V hi -^- + A^B + V s (r)-Ej\^ s (r) 



. (40) 



Here, Vr(r) and Vg(r) are the interaction potentials of atoms with internal 
state |T) and |S), respectively, and A/iB is their difference in Zeeman energy 
due to the interaction with the magnetic field B, with Afj, the difference in 
magnetic moment. In agreement with the above remarks, |T) is referred to 
as the triplet channel, whereas |S) is referred to as the singlet channel. The 
potentials Vr(r) and Vs(r) are the triplet and singlet interaction potentials, 
respectively. 

As a specific example, we use for both interaction potentials again square well 
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potentials, 



Vr, s (r) 



-Vr,s if r < R 
if r > R 



(41) 



where Vr,s > 0. For convenience we have taken the range the same for both 
potentials. Furthermore, we assume that the potentials are such that V T < V s 
and that Vs is just deep enough such that it contains exactly one bound state. 
Finally, we assume that < Vhf <C Vr, Vs, AfiB. The potentials are shown in 
Fig. 5. 

To discuss the scattering properties of the atoms, we have to diagonalize the 
hamiltonian for r > R, in order to determine the incoming channels, which are 
superpositions of the triplet and singlet states |T) and |S). Since the kinetic 
energy operator is diagonal in the internal space of the atoms, we have to find 
the eigenvalues of the hamiltonian 



! V hf N 



\V h{ AfiB J 
These are given by 



(42) 



(43) 



The hamiltonian ff > is diagonalized by the matrix 



Q(0) 



cos 9 sin 9 
- sin 9 cos 9 



(44) 



according to 



Q(9 > )H > Q- 1 (9 > ) 



e> 
e> 



(45) 



which determines tan9 > = — IV^j A\iB . We define now the hyperfine states 
\\\) and III) according to 
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= Q(0>) 




JU>J 




I IS) J 



(46) 



which asymptotically represent the scattering channels. In this basis the Schrodinger 
equation for all r reads 



V 




(47) 



where the energy is measured with respect to e> and we have defined the 
potentials according to 



V n (r)V n (r)\ = /V T (r) 

K V n (r) V u (r) ) { V s (r)J 



(48) 



Since all these potentials vanish for r > R we can study scattering of atoms in 
the states | ft) and | J, J,). Because the hyperfine interaction \4f is small we have 
that ~ A/J.B and e> ~ 0. Moreover, for the experiments with magnetically- 
trapped gases we always have that AfiB ^> k^T where /cb is Boltzmann's 
constant and is T the temperature. This means that in a realistic atomic gas, 
in which the states | ft) and | J, J,) are available, there are in equilibrium almost 
no atoms that scatter via the latter state. Because of this, the effects of the 
interactions of the atoms will be determined by the scattering amplitude in the 
state | If). If two atoms scatter in this channel with energy E ~ ksT AfiB 
they cannot come out in the other channel because of energy conservation. 
Therefore, the indices ft refers to an open channel, whereas J, J, is associated 
with a closed channel. The situation is further clarified in Fig. 5. 

To calculate the s-wave scattering length in the open channel we have to solve 
the Schrodinger equation. In the region r > R the solution is of the from 




Ce ikr + De~ ikr 
Fe~ Kr 



(49) 



where k = y'm(e^ — e>)/h 2 — k 2 and, because we have used the same notation 
as in Eq. (12), the s-wave phase shift is again determined by Eq. (13). In the 
region r < R the solutions are of the form 
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V T (r) 



Fig. 5. Feshbach resonance in a two-channel system with square- well interaction 
potentials. The triplet potential Vr(r) is indicated by the thick dashed line. The 
singlet potential that contains the bound state responsible for the Feshbach res- 
onance is indicated by the thin dashed line. Due to the Zeeman interaction with 
the magnetic field, the energy difference between the singlet and triplet is equal to 
AfiB. The interactions in the open and closed hyperfine channels are indicated by 
Vj^(r) and Vj^(r), respectively. 



(50) 





-( 








where 



k< = ym(e> - e<)/h 2 + k 2 ; 

k< = ^fm(e> - e<)/h 2 + k 2 , (51) 



and 



4 = — ^ Vs Ty(V s -V T -A^B) 2 + (2V M ) 2 . (52) 



are the eigenvalues of the matrix 
/ 



V T V m 
\ V M AfxB-Vs 



(53) 



In order to determine the phase shift we have to join the solution for r < R 
and r > R smoothly. This is done most easily by transforming to the singlet- 
triplet basis {|T), |S)} since this basis is independent of r. Demanding the 
solution to be continuously differentiable leads to the equations 
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4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 4.9 

A uB/(n 2 /mR 2 ) 



Fig. 6. Scattering length for two coupled square-well potentials as a function of 
A[iB. The depth of the triplet and singlet channel potentials is Vt = h 2 /mR 2 and 
Vs = Wh 2 /mR 2 , respectively. The hyperfine coupling is Vhf = 0.1ft 2 /mR 2 . The 
dotted line shows the background scattering length a^ g . 



Q-\e<) 



u 



TT 



u 



11 



(R) 
(R) 



Q-\o>) 



d_ 

dr 



Q 



-i 





(54) 



--R 



where taxi6 < = 2Vhf/(Vs — Vt — AfiB). These four equations determine the 
coefficients A, B, C, D and F up to a normalization factor, and therefore also 
the phase shift and the scattering length. Although it is possible to find an 
analytical expression for the scattering length as a function of the magnetic 
field, the resulting expression is rather formidable and is omitted here. The 
result for the scattering length is shown in Fig. 6, for Vs = 10ft 2 /mi? 2 , Vp = 
h 2 /mR 2 and Vhf = 0.1ft 2 /mR 2 , as a function of AfiB. The resonant behaviour 
is due to the bound state of the singlet potential Vs(r). Indeed, solving the 
equation for the binding energy in Eq. (21) with Vq — ~ we nn d that 
\E m \ ~ 4.62ft 2 /mR 2 , which is approximately the position of the resonance in 
Fig. 6. The difference is due to the fact that the hyperfine interaction leads to 
a shift in the position of the resonance with respect to E m . 

The magnetic-field dependence of the scattering length near a Feshbach reso- 
nance is characterized experimentally by a width AB and position B$ accord- 
ing to 



a(B) = dbg ( 1 - 



AB 
B — Bq 



(55) 



This explicitly shows that the scattering length, and therefore the magnitude 
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-0.02 1 ' ' ' ' 1 

4.615 4.62 4.625 4.63 4.635 4.64 

A uB/(n 2 /mR 2 ) 

Fig. 7. Bound-state energy of the molecular state near a Feshbach resonance for 
two coupled square-well interaction potentials. The solid line and the inset show 
the result for Vhf = 0.1ft 2 /mR 2 . The dashed line corresponds to Vm = 0. The other 
parameters are the same as in Fig. 6. 

of the effective interatomic interaction, may be altered to any value by tuning 
the magnetic field. The off-resonant background scattering length is denoted 
by cibg and is, in our example, approximately equal to the scattering length of 
the triplet potential Vr(r). Using the expression for the scattering length of a 
square well in Eq. (18) for 7 = 1, we find that abg — —0.56/?. Furthermore, 
we have for our example that the position of the resonance is given by B ~ 
4.64h 2 /mAfiR 2 and that the width is equal to AB ~ -0.05ft 2 /mAfiR 2 . 

Next, we calculate the energy of the molecular state for the coupled-channel 
case which is found by solving Eq. (47) for negative energy. In particular, 
we are interested in its dependence on the magnetic field. In the absence 
of the hyperfine coupling between the open and closed channel we simply 
have that e m (B) = E m + AfxB. Here, E m is the energy of the bound state 
responsible for the Feshbach resonance, that is determined by solving the 
single-channel Schodinger equation for the singlet potential. This bound-state 
energy as a function of the magnetic field is shown in Fig. 7 by the dashed line. 
A nonzero hyperfine coupling drastically changes this result. For our example 
the bound-state energy is easily calculated. The result is shown by the solid line 
in Fig. 7 for the same parameters as before. Clearly, close to the resonance the 
dependence of the bound-state energy on the magnetic field is no longer linear, 
as the inset of Fig 7 shows. Instead, it turns out to be quadratic. Moreover, 
the magnetic field B where the bound-state energy is equal to zero is shifted 
with respected to the case where Vhf = 0. It is at this shifted magnetic field 
that the resonance is observed experimentally. Moreover, for magnetic fields 
larger than B there no longer exists a bound state and the molecule now 
decays into two free atoms due to the hyperfine coupling, because its energy 
is above the two-atom continuum threshold. 
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Close to resonance the energy of the molecular state turns out to be related 
to the scattering length by 



UB) = ~^hbW ' (56) 

as in the single-channel case. As we will see in the next sections, the reason for 
this is that close to resonance the effective two-body T-matrix again has a pole 
at the energy in Eq. (56). This important result will be proven analytically in 
Section 4. First, we derive a description of the Feshbach resonance in terms of 
coupled atomic and molecular quantum fields. 



3 Many-body theory for Feshbach-resonant interactions 

In this section we derive the effective quantum field theory that offers a descrip- 
tion of Feshbach-resonant interactions in terms of an atom-molecule hamilto- 
nian. We start from a microscopic atomic hamiltonian that involves atoms 
with two internal states, i.e., we consider a situation with an open and a 
closed channel that are coupled by the exchange interaction. The first step 
is to introduce a quantum field that describes the bound state in the closed 
channel, which is responsible for the Feshbach resonance. This is achieved us- 
ing functional techniques by a so-called Hubbard-Stratonovich transformation 
and is described in detail in Section 3.1. This section is somewhat technical 
and may be omitted in a first reading of this paper. The most important result 
is a bare atom-molecule quantum field theory that is presented in Sec. 3.2. In 
Section 3.3 we subsequently dress the coupling constants of this bare atom- 
molecule theory with ladder diagrams, to arrive at the desired effective quan- 
tum field theory that includes all two-atom physics exactly. The Heisenberg 
equations of motion of this effective field theory are presented in Section 3.4. 

3.1 Bare atom-molecule theory 

Without loss of generality we can consider the simplest situation in which 
a Feshbach resonance arises, i.e., we consider a homogeneous gas of identical 
atoms in a box of volume V. These atoms have two internal states, denoted by 
1 1) and 1 1), that are described by the fields 0j (x, r) and 0j,(x, r), respectively. 
The atoms in these two states interact via the potentials V^(x — x') and 
Vj.j.(x — x'), respectively. The state | j) has an energy AfiB/2 with respect 
to the state | f) due to the Zeeman interaction with the magnetic field B. 
The coupling between the two states, which from the point of view of atomic 
physics is due to the difference in singlet and triplet interactions, is denoted by 
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V|j(x— x'). Putting everything together we write the grand-canonical partition 
function for the gas as a path integral given by 



gr 



(57) 



Since we are dealing with bosons, the integration is over all fields that are 
periodic on the imaginary-time axis ranging from zero to h/3, with h Planck's 
constant and (3 = 1/ k^T the inverse thermal energy. The Euclidian action is 
given by 



h/3 




h{ 


fdx 


v 





s[4>*,<fa,<i>l,<i>i] = 

d d 
*(x, r)H— T (x, r) + 0*(x, r)^— 0|(x, r) 

with the grand-canonical hamiltonian functional given by 



(5* 



H[4>* v (t>h ( t>\, ( t>i]= j dx0*(x,r) 



ft 2 V 2 
' 2m 



+ -|rfx / ^(x',r)\/ TT (x-x / )0 T (x',r)] T (x,r) 
+ /^(x,r) + _^_- / , 

+i / dx'0!(x',r)y u (x - x')<^(x',r)] <^(x,r) 
+ i/rfx/rfx' 



> T (x,Tj(/>^x,Tj 



xF n (x - x , )0 i (x', r)^(x, r) + c.c] , 



(59) 



where \x is the chemical potential of the atoms. Note that this hamiltonian 
functional is the grand-canonical version of the hamiltonian in Eq. (47). The 
indices j and | now refer again to single-particle states, and the two-particle 
hyperfine states are denoted by | ft) an d | 1 1), respectively. The closed- 
channel potential is assumed again to contain the bound state responsible for 
the Feshbach resonance, as illustrated in Fig. 8. 

As a first step towards the introduction of the molecular field that describes the 
center-of-mass motion of this bound state, we introduce the complex pairing 
field A(x, x', r) and rewrite the interaction in the closed channel as a gaussian 
functional integral over this field, given by 
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closed channel 




A |i B 



open channel 



Fig. 8. Illustration of a Feshbach resonance. The upper potential curve corresponds 
to the closed-channel interaction potential Vj^(x — x') that contains the bound state 
responsible for the Feshbach resonance, indicated by the dashed line. The lower 
potential curve corresponds to the open-channel interaction potential Vj-|-(x — x'). 

exp j--^ J dx'0J(x, r)0J(x', r)V n (x - x')0|(x', r)^(x, r) 



propto / d[A*]d[A] exp < 



h/3 



2h 



dr / <ix / rfx' 



A*(x, x / ,r)0 i (x / ,r)0 i (x,r) 



+0J(x', r)0J(x, r)A(x, x', r) - A*(x, x', t)V^(x - x')A(x, x', r) 



• (60) 



This step is known as a Hubbard-Stratonovich transformation [95,96] and 
decouples the interaction in the closed channel. In the BCS-theory of super- 
conductivity this Hubbard-Stratonovich transformation introduces the order 
parameter for the Bose-Einstein condensation of Cooper pairs into the theory. 
This order parameter is the macroscopic wave function of the condensate of 
Cooper pairs. We shall see below that, in our case, the role of the Cooper pair 
is played by the diatomic molecular state that is responsible for the Feshbach 



resonance. 



The functional integral over the fields </>*(x, r) and 0^(x, r) has now become 
quadratic and we write this quadratic part as 



h 



h/3 



h(3 



dr dx dr' / dx' 



0I(x,r),^(x,r) 



G u (x,r;x,r) 



0|(x',r') 
0t(x',r') 



(61) 
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where the so-called Nambu-space Green's function for the closed channel obeys 
the Dyson equation 



GTftx, r; x', t') = G -|, (x, r; x', r') - £ u (x, r; x', r') . 



(62) 



The noninteracting Nambu-space Green's function is given by 



G ,u(x,t;x\t') 
where 



G o,li( x > r ; x, > r ') 



G - a \(x',r';x,r) 



(63) 



9 h 2 V 2 AfiB 



G 0)U (x, r; x\ t') = -H5(r - t')5(x - x')(64) 



is the single-particle noninteracting Green's function. The self-energy is purely 
off-diagonal in Nambu space and reads 



nS u (x,r;x / ,r / ) = 5(r-r / ) 
where 



k(x, x',t) 

k*(x,x',t) 



K 



(x, x', r) = A(x, x', r) + % (x - x')0 T (x, r)0 T (x', r) . 



(65) 



(66) 



Note that a variation of the action with respect to the pairing field shows that 
(A(x,x',r)> = (y u (x-x')0i(x)^(x')> , (67) 



which relates the auxiliary pairing field to the wave function of two atoms 
in the closed channel. Roughly speaking, to introduce the field that describes 
a pair of atoms in the closed-channel bound state we have to consider only 
contributions from this bound state to the pairing field. Close to resonance it is 
this contribution that dominates. Note that the average of the pairing field in 
Eq. (67) indeed shows that the pairing field is similar to the macroscopic wave 
function of the Cooper-pair condensate. However, in this case we are interested 
in the phase (A) = and therefore need to consider also fluctuations. 

Since the integration over the fields (x, r) and <fi ^ (x, r) involves now a gaus- 
sian integral, it is easily performed. This results in an effective action for the 
pairing field and the atomic fields that describes the open channel, given by 
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S cS i<P*, T , A*, A] = Jdr Jdx j 0*(x, r)h— T (x, r) 



n 2 v 2 

' 2m 



+ dx'0*(x', r)y TT (x - x')0 T (x', r) 



+ <^(x,r) 

~\l dT / rfx ' [ A *( X ' x/ ' T ) V ii 1 ( X - x/ ) A ( X > X » 



h 



T (x,r) 



> H-g:, 1 ) 



(68) 



Because we are interested in the bare atom-molecule coupling we expand the 
effective action up to quadratic order in the fields A*(x, x', r) and A(x, x', r). 
Considering higher orders would lead to atom-molecule and molecule-molecule 
interaction terms that will be neglected here, since in our applications we 
always deal with a small density of molecules relative to the atomic density. 



Hence, we expand the effective action by making use of 

OO 1 

Tr[ln(-G7/)] = TrM-Goi)] - £ -Tr[(G , u S u ) m ] 



(69) 



m=l 



This results for the part of the effective action that is quadratic in A*(x, x', r) 
and A(x, x', r) in 



1 



hp 



HP 



S[A*,A] = --Jdr Jdx Jdx' Jdr' Jdy J dy' 
o o 
x A*(x, x', t)HG?(x, x', r; y, y', r') A(y, y', r') , 

where the Green's function of the pairing field obeys the equation 



(70) 



G A (x, x\ r; y, y', r') = hV u (x - x')5(x - y)5(x' - y')5(r - r') 

hp 

\ J dr" JdzJ dz' [F u (x - x')G 0lU (x, r; z, t")G , u (x', r; z', r") 
o 

xG A (z,z',r";y,y',r')] . 



(71) 



From this equation we observe that the propagator of the pairing field is related 
to the many-body T-matrix in the closed channel. More precisely, introducing 
the Fourier transform of the propagator to relative and center-of-mass mo- 
menta and Matsubara frequencies Vt n = 2irn/ h(3, denoted by G A (k, k', K, iVt n ), 
we have that 



G A (k, k', K, iVL n ) = hTfl B (k, k', K, itel n - AfiB + 2/i) , 



(72) 
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where the many-body T-matrix in the closed channel obeys the equation 



T^ B (k,k',K,z) = \/ u (k-k') 



1 + N (e K / 2+ k»-/i + ^) + N ( eK/2 _ k „- /i +^) 



V k „ Z — CK/2+k" — CK/2-k" 

Ml 

U 



7™ B (k",k',K,z) . (73) 



with N(x) = [e^ x — the Bose distribution function. Here, V^(k) = 
J rfxVjj(x)e ik ' x denotes the Fourier transform of the atomic interaction po- 
tential. This equation describes the scattering of a pair of atoms from relative 
momentum k' to relative momentum k at energy z. Due to the fact that the 
scattering takes places in a medium the many-body T-matrix also depends 
on the center-of-mass momentum K, contrary to the two-body T-matrix in- 
troduced in the previous section, which describes scattering in vacuum. The 
kinetic energy of a single atom is equal to = h 2 k. 2 /2m. The factor that 
involves the Bose-Einstein distribution function arises because the probability 
of a process where a boson scatters into a state that is already occupied by 
Ni bosons is proportional to 1 + i\q. The reverse process is only proportional 
to Ni. This explains the factor 



1 + N 1 + N 2 = (1 + N 1 )(1 + N 2 )-N 1 N 2 , (74) 
in the equation for the many-body T-matrix [97]. 

The many-body T-matrix is discussed in more detail in the next section when 
we calculate the renormalization of the interatomic interactions. For now we 
only need to realize that, for the conditions of interest to us, we are always 
in the situation where we are allowed to neglect the many-body effects in 
Eq. (73) because the Zeeman energy Afj,B/2 strongly suppresses the Bose oc- 
cupation numbers for atoms in the closed channel. This is certainly true for 
the experimental applications of interest because in the current experiments 
with magnetically-trapped ultracold gases the Zeeman splitting of the mag- 
netic trap is much larger than the thermal energy. This reduces the many-body 
T-matrix equation to the Lippmann-Schwinger equation in Eq. (29) for the 
two-body T-matrix in the closed channel T, 2 , B (k, \s!,z — ck/2), which, in its 
basis-independent operator formulation, reads 



1 

z~H 



fff(z) = V U + V U ^-fffiz) , (75) 



with H = p 2 /m. As we have seen previously, this equation is formally solved 
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by 



rp2B. 

J u 1 



z-H 



-V\ 



U ' 



(76) 



II 



with = H + Vj^. From the previous section we know that the two-body 
T-matrix has poles at the bound states of the closed-channel potential. We 
assume that we are close to resonance and hence that one of these bound 
states dominates. Therefore, we approximate the two-body T-matrix by 

Tff (,) ^ u %^/ u , (77) 



where the properly normalized and symmetrized bound-state wave function 
X m (x) = (x|% m ) obeys the Schrodinger equation 



h 2 v 2 



m 



+ Vu(*) 



Xm(x) = ^mXm(x) . 



(78) 



It should be noted that this wave function does not correspond to the dressed, 
or true, molecular state which is an eigenstate of the coupled-channels hamilto- 
nian and determined by Eq. (47). Rather, it corresponds to the bare molecular 
wave function. The coupling V^(x — x') of this bare state with the continuum 
renormalizes it such that it contains also a component in the open channel. 
Moreover, as we have already seen in the previous section, this coupling also 
affects the energy of this bound state. Both effects are important near the 
resonance and are discussed in detail later on. 

We are now in the position to derive the quadratic action for the quantum field 
that describes the bare molecule. To do this, we consider first the case that the 
exchange interaction V^(x — x') is absent. Within the above approximations, 
the two-point function for the pairing field is given by 

<A< k ,K,«„)A>( k <,K,a,)> = .(79) 



We introduce the field m (x, r), that describes the bound state in the closed 
channel, i.e, the bare molecule, by considering configurations of the pairing 
field such that 

A(x,x',r) = v / 2ni(x-x , )x m (x-x')0 m ((x + x')/2,r) . (80) 
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Using this we have that 

<«k, a„, C (K, a.)) = _, mn + lK/2+ h Em + A,B-2, • (81 > 

which shows that the quadratic action for the bare molecular field is, in posi- 
tion representation, given by 

h/3 

^[0m> 0m] = Jdr J dx (f)* m (x, t) 



m (x,r). (82) 

In the absence of the coupling of the bare molecular field to the atoms, the 
dispersion relation of the bare molecules is given by 

twk(B) = e k /2 + E m + AfiB . (83) 



x 



d h 2 v 2 



As expected, the binding energy of the bare molecule is equal to e m (B) = 
E m + AfiB. The momentum dependence of the dispersion is due to the kinetic 
energy of the molecule. 

To derive the coupling of this bare molecular field to the fields </>* (x, r) and 
0j(x, r) it is convenient to start from the effective action in Eq. (68) and to 
consider again only terms that are quadratic in the self-energy. Integrating 
out the pairing fields leads to an interaction term in the action for the field 
describing the open channel, given by 



hf3 hf3 
l -Jdr jdx Jdx'Jdr'Jdy Jdy' [v n (x - x')^(x, r)^(x', r 



o 



x Gi 4 i ) (x,x',r;y,y',r')V n (y-y')0 T (y,r')0 T (y / ,r')l , (84) 



where the two-atom four-point Green's function is given diagrammatically in 
Fig. 9. For our purposes it is, for the same reasons as before, sufficient to 
neglect the many-body effects on this propagator and to consider again only 
the contribution that arises from the bound state in the closed channel. This 
gives for the Fourier transform of this Green's function 

G {4) (k k' K Q ) X*n( k )Xm(k') 

G ^ (k ' k ' K ' * ln) ~ ihQ n - e K /2 - AmB - E m + 2/* ' { ^ } 
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Fig. 9. Diagrammatic representation of the two-particle Green's function in the 
closed channel. The solid lines correspond to single-atom propagators. 

where Xm(k) is the Fourier transform of the bound-state wave function. Af- 
ter substitution of this result into Eq. (84) the resulting interaction term is 
decoupled by introducing the field m (x, r) with the quadratic action given 
in Eq. (82). This procedure automatically shows that the bare atom- molecule 
coupling constant is equal to V^(k)x m (k)/-\/2. 



3.2 Bare atom-molecule hamiltonian 



In the previous section we have derived, from a microscopic atomic hamilto- 
nian, a bare atom-molecule theory for the description of a Feshbach resonance. 
It is determined by the action 



hp . 
JdrUdx 



o 



d d 
<j>*(x, r)fr— T (x, r) + </4(x, r)h— m (x, r) 



dr 



'p VP <fm, VmJ f ) 



(86) 



where the bare or microscopic atom-molecule hamiltonian functional is given 
by 



H[(p* v (pp(j) m ,(j)* m ] = 
dtx.(f)* (x, r) 



2m 

h 2 v 2 



^ + \J dx>j(x', r)\/ n (x - x')0 T (x', r) 



T (x,r) 



4m 



+ AfiB + E m -2/2 



(x,r) 



+ jdx |dx'[ 5n (x-x')C((x + x')/2,r)0 T (x',r)0 T (x,r)+c.c] , (87) 

and the bare atom-molecule coupling is given by ^(x) = V|j(x)x m (x)/y / 2, 
where Vjj_(x) is the coupling between the open and closed atomic collision 
channel of the Feshbach problem, that has its origin in the exchange interaction 
of the atoms. Note also that the atom-molecule coupling is proportional to 
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the wave function ^ m (x) for the bound molecular state in the closed channel 
responsible for the Feshbach resonance. 



Physically, the microscopic hamiltonian in Eq. (87) describes bosonic atoms 
in the open channel of the Feshbach problem in terms of the fields 0?(x, r) 
and </>f(x, t). These atoms interact via the interaction potential V^(x — x'). 
Apart from this background interaction, two atoms in the gas can also form a 
molecular bound state in the closed channel with energy E m that is detuned 
by an amount of AfiB from the open channel. This bare molecular state is 
described by the fields 0^(x, r) and m (x, r). The most important input in 
the derivation of Eq. (87) is that the energy difference between the various 
bound states in the closed channel is much larger than the thermal energy, so 
that near resonance only one molecular level is of importance. This condition 
is very well satisfied fo almost all the atomic gases of interest. An exception is 
6 Li, which has two Feshbach resonances relatively close to each other [98,63]. 
The derivation presented in the previous section is easily generalized to this 
situation, by introducing an additional molecular field to account for the sec- 
ond resonance. 



To point out the differences of our approach with work of other authors a 
few remarks are in order. First of all, our starting point was the microscopic 
two-channel atomic hamiltonian in Eq. (59), from which we derived the mi- 
croscopic atom- molecule hamiltonian in Eq. (87). As we started with the full 
interatomic interaction potentials, the atom-molecule coupling constant and 
atom-atom interaction have momentum dependence which cut off the momen- 
tum integrals encountered in perturbation theory. Because of this, no ultravio- 
let divergencies are encountered at any order of the perturbation theory, as we 
will see in the next section. This contrasts with the model used by Kokkelmans 
and Holland [79], and Mackie et al. [78], who use a phenomenological atom- 
molecule hamiltonian with delta-function interactions and therefore need a 
renormalization procedure to subtract the ultraviolet divergencies. 



In an application of the above microscopic atom-molecule hamiltonian to 
realistic atomic gases we have to do perturbation theory in the interaction 
V||(x — x') and the coupling g^(x — x'). Since the interatomic interaction is 
strong, this perturbation theory requires an infinite number of terms. Progress 
is made by realizing that the atomic and molecular densities of interest are 
so low that we only need to include two-atom processes. This is achieved by 
summing all ladder diagrams as explained in detail in the next section. 
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3.3 Ladder summations 



From the bare or microscopic atom-molecule theory derived in the previous 
section we now intend to derive an effective quantum field theory that con- 
tains the two-atom physics exactly. This is most conveniently achieved by 
renormalization of the coupling constants. Moreover, the molecules acquire 
a self-energy. Both calculations are done within the framework of perturba- 
tion theory to bring out the physics involved most clearly. It is, however, also 
possible to achieve the same goal in a nonperturbative manner by a second 
Hubbard-Stratonovich transformation. 

Because we are dealing with a homogeneous system, it is convenient to perform 
the perturbation theory in momentum space. Therefore, we Fourier transform 
to momentum space, and expand the atomic and molecular fields according 
to 



T (x,r) 



172 a k,> 



Ak-x—iuj n T 



( W/ 2 £ 



(88) 



and 



<Mx, r) = £ 6 k , n e ik — , (89) 

respectively. The even Matsubara frequencies u n = 2im/hf3 account for the 
periodicity of the fields on the imaginary-time axis. With this expansion, the 
grand-canonical partition function of the gas is written as a functional integral 
over the fields ak, n and b^ n and their complex conjugates. It is given by 



Z gr = J d[a*]d[a]d[b*}d[b}exp{--S[a*,a,b*,b}} , (90) 

where the action S[a*,a,b*,b] is the sum of four terms. The first two terms 
describe noninteracting atoms and noninteracting bare molecules, respectively, 
and are given by 

S a [a*, a] = {-iftujn + e k - //) a* Kn a^ n , (91) 

k,n 



and 

Sm[b\ 6] =E H^n + e k/2 + E m + AfiB - 2n) b* Kn b Kn . (92) 

k,n 
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Fig. 10. (a) Ladder diagrams that contribute to the renormalization of the inter- 
atomic interaction, (b) Diagrammatic representation of the Lippmann-Schwinger 
equation for the many-body T-matrix. The solid lines correspond to single-atom 
propagators. The wiggly lines correspond to the interatomic interaction V^. 



The atomic interactions are described by the action 



S int [a*, a] — -— — ^2 ^tt(k — k')aK/2+k,n/2+m a K/2-k,n/2-m 
Zn l JV K,k,k' 

n,m,m' 

—m! i 

(93) 

where Vfj (k) is the Fourier transform of the interatomic interaction potential. 
This Fourier transform vanishes for large momenta due to the nonzero range 
of the interatomic interaction potential. The last term in the action describes 
the process of two atoms forming a molecule and vice versa, and is given by 



S coup [a*,a,b*,b] 



0/2 E<m( k ) 



K.k 



6K il i«K/2+k 1 i!/2+m«K/2-k,n/2-m + C.C. 



X 



(94) 



where <?i-j.(k) is the Fourier transform of the bare atom- molecule coupling 
constant. This coupling constant also vanishes for large momenta since the 
bare molecular wave function has a nonzero extent. 



We first discuss the renormalization of the microscopic atomic interaction 
V||(k), due to nonresonant background collisions between the atoms. The first 
term that contributes to this renormalization is of second order in the interac- 
tion. It is found by expanding the exponential in the path-integral expression 
for the grand-canonical partition function in Eq. (90). To second order in the 
interactions this leads to 



Z gr = J d[a*]d[a] (l - -S int [a*,a\ + ^S^ nt [a*,a] + 



x exp 



h 



S a [a*,a] 



(95) 
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After the decoupling of the eight-point function resulting from the square of 
the action Si nt [a*,a] with the use of Wick's theorem, it gives rise to various 
terms in the perturbation theory which can be depicted by Feynman diagrams 
[96,99]. As mentioned already, we only take into account the ladder Feynman 
diagram. This diagram is given by the second term of the Born series depicted 
in Fig. 10 (a), and corresponds to the expression 



)- ]T Vf T (k - k")G ,a (K/2 + k", iu n/2+m ) 



xG , a (K/2 - k", iu n/2 - m ) V n (k" - k') , (96) 

where 

G , a (k, iw n ) = — ^ , (97) 



is the noninteracting propagator of the atoms. After performing the sum- 
mation over the Matsubara frequencies we find that, to second order, the 
renormalization of the interatomic interactions is given by 



Vtt(k-k') -> Vt T (k-k') 



1 + N (e K / 2 +k" -ft) + N (e K /2-k» -/•*) 



+ - VV^k-k", 

V k „ inw n — e K /2+k" — eK/2-k" + 2/i 

xy n (k"-k'), (98) 

which is finite due to the use of the true interatomic potential. In comparing 
this result with the first two terms of the Born series for scattering in vacuum 
in Eq. (31), we see that the only difference between the two-body result and 
the above result is the factor involving the Bose distributions. This so-called 
statistical factor accounts for the fact that the scattering takes place in a 
medium and is understood as follows. The amplitude for a process where 
an atom scatters from a state with occupation number N\ to a state with 
occupation number N 2 contains a factor N\(l + iV 2 ). The factor N\ simply 
accounts for the number of atoms that can undergo the collision, and may be 
understood from a classical viewpoint as well. However, the additional factor 
(1 + N 2 ) is a result of the Bose statistics of the atoms and is therefore called 
the Bose-enhancement factor. For fermions this factor would correspond to the 
Pauli-blo eking factor (1 — N 2 ), reflecting the fact that a fermion is not allowed 
to scatter into a state that is already occupied by an identical fermion. In 
calculating the Feynman diagram we have to take into account the forward 
and backward scattering processes, which results in the statistical factor in 
Eq. (98). 
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MB 



Fig. 11. Renormalization of the atom- molecule coupling constant by interatomic 
interactions. The solid lines correspond to single-atom propagators. The wiggly lines 
corresponds to the interatomic interaction V^. 

Continuing the expansion in Eq. (95) and taking into account only the ladder 
diagrams leads to a geometric series, which is summed by introducing the 
many-body T-matrix in the open channel. It is given by 



T T 7 B (k,k',K,z) = \/ TT (k-k') 



+^E^tt( k - k ") 



1 + N (e K /2+k" - A*) + N (e K /2-k" ~A*) 



V k „ Z — CK/2+k" — e K/2-k" 

iM] 

TT 



xT™ B (k",k',K,z) . (99) 



Its diagrammatic representation is given in Fig. 10 (b). For the moment we 
neglect the many-body effects on the scattering atoms and put the Bose- 
distribution functions equal to zero. This assumption is valid at temperatures 
far below the critical temperature [100]. This reduces the many-body T-matrix 
to the two-body T-matrix T?P(k, k', z — ck/2). For the low temperatures of 
interest to us here, we are allowed to take the external momenta equal to 
zero. For small energies we find, using the result in Eq. (38), that the effective 
interaction between the atoms reduces to 



T T 2 f(0, 0, ihuj n - e K /2 + 2/x) 



■TT 



m 



x 



1 



i _ / -m(ifrw n -E K /2+2ji) _ ab g r bg m(^ n -e K /2+2^) 

1 "bgV fr2 2h 2 



(100) 



Here ab g and r bg are the scattering length and the effective range of the open- 
channel potential V^{x), respectively. Although these could in principle be 
calculated with the precise knowledge of this potential, it is much easier to 
take them from experiment. For example, the magnitude of the scattering 
length can be determined by thermalization-rate measurements [4]. The effec- 
tive range is determined by comparing the result of calculations with experi- 
mental data. We will encounter an explicit example of this in Section 6. 

The next step is the renormalization of the microscopic atom-molecule cou- 
pling constant. Using the same perturbative techniques as before, we find that 
the effective atom-molecule coupling is given in terms of the bare coupling by 
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g MB (k, K, z) = g n (k) + 1 £ T^ IB (k, k', K, z) 



x 



1 + iV(eK/2+k' - A-0 + N(e K / 2 -k> - jj) 



Z — e K/2+k' — e K/2-k' 



(101) 



and is illustrated diagrammatically in Fig 11. Neglecting again many-body 
effects, the coupling constant becomes g 2B (k, z — ck/2) with 



/ B (k^)=^ n (k)+i5:T T 2 T B (k,k'^) 



z — 2ek 



■9nM ■ 



(102) 



From the above equation we infer that the energy dependence of this coupling 
constant is the same as that of the two-body T-matrix. This result is easily 
understood by noting that for a contact potential Vfj(k) = V and we simply 
have that g 2B = g^T^ /Vq. Hence we have for the effective atom-molecule 
coupling 



g 2B (0,thw n -e K /2 + 2ii) 



9 



i _ / -m(ihhj n -e K /2+2jj!) 

[i - a bg V 



a bg r bg m(ifbuj n -e- K /2+2Li) 
2h 2 



(103) 



where g is the effective atom-molecule coupling constant at zero energy. The 
latter is also taken from experiment. We come back to this point in Section 4.1 
where we discuss the two-atom properties of our effective many-body theory. 

Finally, we have to take into account also the ladder diagrams of the resonant 
part of the interaction. This is achieved by including the self-energy of the 
molecules. It is in first instance given by the expression 



n MB (K,,) = ^X>n(k) 



xg M *(k,K,z) , 



1 + iV(e K /2+k - //) + iV(e K / 2 -k - aO 



Z — CK/2+k — ^K/2-k 



(104) 



and shown diagrammatically in Fig. 12. We neglect again many-body effects 
which reduces the self-energy in Eq.(104) to I1 2B (2; — £k/2) with 



n^(z) = ( Xm \V n G n (z)V n \ Xin ) , 



(105) 



where the propagator G^(z) is given by 
1 



G n (z) 



z-H 



(106) 



TT 
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Fig. 12. Molecular self-energy. The solid lines correspond to single-atom propaga- 
tors. The wiggly lines corresponds to the interatomic interaction V^. 

with the hamiltonian 



= h V\ 



m 



TT 



TT 



(107) 



We insert in Eq. (105) a complete set of bound states \ip K ) with energies E K 
and scattering states that obey the equation in Eq. (25). This reduces 

the self-energy to 



n 2B (^) = EI(xm|Vul<M 



z - E K 



(108) 



+ 



/ 



dk 



l<Xm|Vul4 +) >| 2 - 



z-2e k ' 



where we replaced the sum over the momenta k by an integral. Using Eq. (102) 
and the equation for the scattering states we have that 



9 2B (kM) = ^(x m \v n \^ +) ) ■ 



(109) 



Neglecting the energy dependence due to the contribution of the bound states 
since their binding energies are always large compared to the thermal en- 
ergy, we have, using the result for to the atom-molecule coupling constant in 
Eq. (103), the intermediate result 



U m (z) = 2 



l/k s 2B (o,24) 2 



(2k) 



z - 2e k 



(110) 



The remaining momentum integral yields the final and for our purposes very 
important result 



hY™{z) = n 2B (^) - n 2B (o) ee rt 2B (z) + (AmB + E m ) 



t2B, 



t2B 



2B, 
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1 + a bg (a b g-r b g) + 



-i 



(111) 



where we have denoted the energy-independent shift n 2B (0) in such a manner 
that the position of the resonance in the magnetic field is precisely at the 
experimentally observed magnetic-field value Bq. This shift is also shown in 
the results of the calculation of the bound-state energy of the coupled square 
wells in Fig. (7). 



3.4 Effective atom-molecule theory 



Putting the results from the previous section together, we find that the atom- 
molecule system is described by the effective action 



S cS [a*, a,b*,b] = J2 (~^n + - aO o*^ n a Kr, 



k,n 



+\jr^ E T S (i^n " W2 + 2/1) 



K,k,k' 
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Xa K/2+k,n/2+m a K/2-k,n/2-m a K/2+k',n/2+m' a K/2-k',n/2-m' 

+ £ [-ihuu n + e k /2 + 5(B) - 2fi + hZ 2B (thuj n - e k /2 + 2/i)] b* Kn b Kn 



k,n 
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— ]T 9 2B (i?^n - e K /2 + 2/i) 



x 



^K,n a K/2+k,n/2+mflK/2-k,ri/2-m + CC. , (H2) 



K,k 

n,m 



where 5(B) = A/j,(B — Bq) is the so-called detuning. From now on we use 
the notation T^(z) = T T 2 T B (0, 0, z), and g 2B (z) = g 2B (0,z). Since these cou- 
pling constants are the result of summing all ladder diagrams, these diagrams 
should not be taken into account again. In the next section we discuss how 
the coupling constants are determined from experiment. 

To consider also the real-time dynamics of the system we derive the Heisenberg 
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equations of motion for the field operators ^ a (x, t) and V> m (x, t), that anni- 
hilate an atom and a molecule at position x and time t, respectively. Their 
hermitian conjugates are the creation operators. To determine the Heisen- 
berg equations of motion for these field operators, we first have to perform an 
analytic continuation from the Matsubara frequencies to real frequencies. To 
ensure that the physical quantities and equations of motion are causal, this 
has to be done by a so-called Wick rotation. This amounts to the replacement 
of the Matsubara frequencies by a frequency with an infinitesimally small and 
positive imaginary part 

iuj n — > uj + . (113) 



This leads to a subtlety involving the analytic continuation of the square root 
of the energy in the various expressions. Due to the branch cut in the square 
root we have that 

\J -ihco n —> yf — (huj + ) = -iVfiw . (114) 



The last expression on the right-hand side of this equation is valid for Hco on 
the entire real axis. 

To obtain the equation of motion in position and time representation, we have 
to Fourier transform back from momentum and frequency space. This amounts 
to the replacement 

d h 2 v 2 

hw-e K /2^ih- + -—. (115) 
at Am 

Note that this combination of time and spatial derivatives is required due to 
the Galilean invariance of the theory. 

For simplicity we assume that we are so close to resonance that we are allowed 
to neglect the energy dependence of the effective atomic interactions and the 
effective atom-molecule coupling. Moreover, for notational convenience we take 
only the leading-order energy dependence of the molecular self-energy into 
account. Higher orders are straightforwardly included but lead to somewhat 
complicated notations in the position and time representation. The leading- 
order energy dependence of the self-energy is, after the Wick rotation to real 
energies, given by 

3 /2 

h^\E) ~ -f^iyfE. (116) 

The additional superscript indicates that we are dealing with the retarded 
self-energy, i.e., the self-energy evaluated at the physically relevant energies 



39 



E + so that KE^(E) = hT?^(E + ). Note that for positive energy E this result 
is in agreement with the Wigner-threshold law. This law gives the rate for 
a state with well-defined positive energy to decay into a three-dimensional 
continuum. 

Within the above approximations, the Heisenberg equations of motion for the 
coupled atom-molecule model read 



ih 



h 2 V 2 Ana h „h 2 



2m 
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dt 

+2#l(x,t)^ m (x,t) , 
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$(x,tty a (x,i) 
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3/2 . L d h 2 V 2 



27chn ik dt + Am 



^a(x,t) 



V> m (x, t) + #a( x , > 



(117) 



where we have also allowed for a time-dependent detuning. In the next sec- 
tion we show that these equations correctly reproduce the Feshbach-resonant 
scattering amplitude and the binding energy of the molecule. Moreover, we 
apply the effective theory derived in this section to study many-body effects 
on this binding energy, above the critical temperature for Bose-Einstein con- 
densation. 



4 Normal state 



In this section we discuss the properties of the gas in the normal state. In the 
first section, we consider the two-atom properties of our many-body theory. 
Hereafter, we discuss the equilibrium properties that follow from our theory. In 
the last section, we investigate many-body effects on the energy of the molec- 
ular state, above the critical temperature for Bose-Einstein condensation. 



4-1 Two-atom properties of the many-body theory 



In this section we show that our effective field theory correctly contains the 
two-atom physics of a Feshbach resonance. First, we show that the correct 
Feshbach-resonant atomic scattering length is obtained after the elimination 
of the molecular field. Second, we calculate the bound-state energy and show 
that it has the correct threshold behaviour near the resonance. To get more 
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insight in the nature of the molecular state near resonance, we also investigate 
the molecular density of states. 



4-1.1 Scattering properties 

To calculate the effective interatomic scattering length, we have to eliminate 
the molecular field from the Heisenberg equations of motion in Eq. (117). Since 
the scattering length is related to the scattering amplitude at zero energy and 
zero momentum, we are allowed to put all the time and spatial derivatives 
in the equation of motion for the molecular field operator equal to zero. This 
equation is now easily solved, which leads to 

Vv(x,t) = --n-^(x,t). ( 118 ) 



Substitution of this result into the equation for the atomic field operator leads 
for the interaction terms to 



*)&(x, t)^ a (x, t) + 2#t( x , i)^ m (x, t) = 

m 

~ W)) ^ i(x ' t)4(x ' t)4(x> t] • (119) 

From this result we observe that we have to take the ren ormalized atom- 
molecule coupling constant at zero energy equal to g — h^jm^ABAp/m, 
so that we have 

47ra bg ft 2 2g 2 Ana(B)h 2 

m 5(B) ~ m ' 1 UJ 



where we recall that the scattering length near a Feshbach resonance is given 
by 

a(B) = a bg - B ^* B ^ j = Q bg + Qres(-B)- (121) 



Since both the width AB and the background scattering length Ob g are known 
experimentally, the knowledge of the difference in magnetic moment between 
the open and the closed channel A/i completely determines the renormalized 
coupling constant g. Since the open and the closed channel usually correspond 
to the triplet and singlet potential, respectively, we always have that |A//| ~ 
2/jL-q, with /xb the Bohr magneton. More precise values of the difference in 
magnetic moments are obtained from coupled-channels calculations using the 
interatomic interaction potentials [10,14,79,101]. 
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From the above analysis we see that the correct Feshbach-resonant scattering 
length of the atoms is contained in our theory exactly. Next, we show that our 
effective theory also contains the correct bound-state energy. 



4-1.2 Bound-state energy 

The energy of the molecular state is determined by the poles of the retarded 
molecular propagator GW(k,w). It is given by 

G^(k,io) = r- . (122) 

huj+ - e k /2 - 5(B) - KE^ifku - e k /2) 



For positive detuning 5(B) there only exists a pole with a nonzero and negative 
imaginary part. This is in agreement with the fact that the molecule decays 
when its energy is above the two-atom continuum threshold. The imaginary 
part of the energy is related to the lifetime of the molecular state. For negative 
detuning the molecular propagator has a real and negative pole correspond- 
ing to the bound-state energy. More precisely, in this case the poles of the 
molecular propagator are given by hio = e m (B) + ek/2, where the bound-state 
energy is determined by solving for E in the equation 

E - 5(B) - HX£\E) = 0. (123) 



In general, this equation cannot be solved analytically but is easily solved 
numerically, and in Section 6 we discuss its numerical solution for the pa- 
rameters of 85 Rb. Close to resonance, however, we are allowed to neglect the 
effective range of the interactions. This reduces the retarded self-energy of the 
molecules to 



t\a hg \ 



Moreover, the bound-state energy is small in this regime and we are allowed to 
neglect the linear terms in the energy with respect to the square-root terms. 
This reduces the equation for the bound-state energy in Eq. (123) to 



,2 



3 / 2 iy/E 



2 ^ 3 l-i|a bg | 



5(B) . (125) 



This equation is easily solved analytically, and yields the result 

UB) = , (126) 

m[a(B)\ 2 
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which analytically proves the numerical result in Eq. (56). This numerical re- 
sult was obtained for the specific case of two coupled attractive square wells. 
The above analytic proof, which does not depend on the details of the poten- 
tial, shows that the result is general. 

The same result is found by noting that after the elimination of the molecular 
field the effective on-shell T-matrix for the atoms in the open channel is given 
by 

T m {E + ) = T b 2 g B (E+) + ^ 2B (£+)| 2 GL +) y^m^E) . (127) 



Close to resonance this expression reduces to 



T (E) ~ 



4:ira rcs (B)h 2 



m 



[l + ia ies (B)J^ 



'128) 



The pole of this T-matrix, which gives the bound-state energy, is indeed equal 
to the result in Eq. (126) close to resonance. 



4-1.3 Molecular density of states 



The molecular density of states is obtained by taking the imaginary part of 
the retarded molecular propagator [99], i.e., 



p m QL,u>) = ~1m [GL +) (k,w) 



(129) 



For simplicity, we discuss here only the situation that we are close to resonance, 
and therefore approximate the retarded molecular self-energy by the square- 
root term resulting from Wigner's threshold law as given in Eq. (116). The 
extension to situations further of resonance are straightforward. 

For the case of negative detuning, the molecular density of states is shown 
by the solid line in Fig. 13 and has two contributions. One arising from the 
pole at the bound-state energy and the second from the two-atom continuum. 
Within the above approximation, it is given by 



Pm(k,w) = Z(B)5(huj - e k /2 - e m (B)) 
+-9(huj - e k /2) 

7T 



X [fkj- e k /2 - 5(B)} 2 + (3 4 m 3 /47r 2 fr 6 ) (Hiv - e k /2) ' { } 
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Fig. 13. Molecular density of states. The solid line shows the density of states for 
negative detuning. Since there is a true bound state in this case there is a pole in 
the density of states. For positive detuning the density of states is approximately a 
Lorentzian as shown by the dashed line. 

with Z(B) the so-called wave-function renormalization factor 



Z(B) 



1 - 



1 + 



du 



n -1 



g 2 m 3/2 
47rfry|e m (5)|_ 



-i -1 



(131) 



This factor goes to zero as we approach the resonance and it becomes equal 
to one far off resonance. Physically, this is understood as follows. Far off res- 
onance the bound state of the coupled-channels hamiltonian in Eq. (47), i.e., 
the dressed molecule, is almost equal to the bound state of the closed-channel 
potential and has zero amplitude in the open channel. This corresponds to 
the situation where Z(B) ~ 1. As the resonance is approached, the dressed 
molecule contains only with an amplitude ^Jz(B) the closed-channel bound 
state, i.e., the bare molecule. Accordingly, the contribution of the open channel 
becomes larger and gives rise to the threshold behaviour of the bound-state 
energy in Eq. (126). Of course, the square of the wave function of the dressed 
molecule is normalized to one. This is expressed by the sum rule for the molec- 
ular density of states, 

J d(fku)p m (k,Lu) = 1 . (132) 



In detail, the dressed molecular state with zero momentum is given by 
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|Xm| dressed) = sJz{B)bl\V>) + ]T C k 4«-k|0) • (133) 

k 

Here, the second-quantized operator b\y creates a molecule with zero momen- 
tum. It acts on the vacuum state |0). The bare molecular state is therefore 
given by \xm) = ^o|0). The operator creates an atom with momentum hk 
and hence the coefficient denotes the amplitude of the dressed molecular 
state to be in the open channel of the Feshbach problem. 

To gain more insight in the nature of the dressed molecular state we cal- 
culate the coefficients Ck in perturbation theory. Neglecting the off-resonant 
background interactions and the energy dependence of the atom-molecule cou- 
pling constant, the hamiltonian appropriate for our purposes is, in terms of 
the above operators, given by 

H = H am + H coup , (134) 



with 



bib 



k°k , 



(135) 



and 



H C oup 



Y [&K«K/2+kaK/2-k + h.C. 
^ K,k 



(136) 



The zeroth-order state around which we perturb is the bare molecular state 
\Xm) with energy 5(B). In first order in g the dressed molecular state is given 
by 



|Xm; dressed) = yJz(B)& o \0) + 



8(B) - H a 



-H coup bi\0) 



= ^t |0> + _^ ? _, 



4«-klo) , 



(137) 



where Z(B) = 1 — 0(g 2 ). This result shows that, to first order in g, the 
coefficients are given by 



9 



1 



VS(B)-2e k 



(138) 



We now calculate the wave-function renormalization factor Z(B) in a different 
manner by demanding that the dressed molecular wave function is properly 
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normalized, i.e., 

(x m ; dressed | x m ; dressed) = 1 . 



(139) 



This leads to 



V V [5(B) - 2e k ]< 



= Z{B) 



du 



(140) 



The factor of two corresponds to the two contributions arising from the matrix 
element (0|aka_k<z k '<zL k /|0). From this result we find that the wave-function 
renormalization factor is given by 



Z(B) = 1 + 



d^+\5(B)) 
dui 



1 - 



d^HMB)) 

du 



-i -l 



(141) 



in agreement with the result in Eq. (131) to second order in g. 

Note that the total number of atoms in the dressed molecular state should be 
equal to two. The number of atoms is given by 

N = 2 X^lA)drcssed + Xl( S k«k) dressed , (142) 
k k 



where (• • -)dressed = (x m ; dressed) • • • \xm', dressed). For the number of atoms 
with momentum hit we have that 



a k a k )dressed - ~y — _ ^ , 



from which, with the use of Eq. (140), we find that 

£(4«k)dressed = 2 - 2Z{B) . (144) 
k 



Using (bib k ) 

dressed 

o we have indeed that the total number of atoms 

iV = 2, as required. 

If the magnetic field varies not too rapidly, we are allowed to make an adiabatic 
approximation to the Heisenberg equation of motion for the bare molecular 
field operator in Eq. (117). This amounts to introducing a molecular field 
V'mC*-) t na t annihilates a dressed molecule, i.e., a molecule with internal state 
given by Eq. (133). This is achieved as follows. In frequency and momentum 
space the action for the bare molecular field is given by 
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S[(j)* m , 0m] = J 7^7 J2 ^m( k > U ) ~ C k/ 2 



doj 

W k 

5(B) - fiEW(fiw - e k /2)l m (k,cu) . 



(145) 



Next, we expand this action around the pole of the propagator e m (B). To 
linear order, this yields the result 



dui 



E 



(2?r) V Jz(B) 



" JX ^ [hu; - e k /2 - e m (B) 



i m (k,cj) 



(146) 



From this equation we see that the field that describes the dressed molecule is 



given by <p' m = <p m / y Z(B). This leads to the following action for the dressed 
molecular field in position and time representation 



s[<CM = j dt j d * C(x,t) 



'in— + — 

dt Am 



e m (B) 



C(x,t). (147) 



More importantly, the terms that describe the coupling between the atoms 



and the molecules are multiplied by a factor yZ(B). In detail, the coupled 
Heisenberg equations of motion for the atomic and dressed molecular field 
operators are given by [41] 



<9^ a (x,t) 



,, ^ m (x,t) 

%h dt 



^ + ^t (x , t)4(x , t) 



4(x,t) 



+2g^t)^,t)^ m (^t) , 

h 2 v 2 



Am 



+ e m (t) 



^,t)+gJZ(t)^,t) , 



(148) 



where Z(t) = Z(B(t)), and e m (t) = e m (B(t)). In the derivation of the above 
coupled equations we have assumed that we are allowed to make an adiabatic 
approximation for the renormalization factor Z{B) and that we can evaluate 
it at every time at the magnetic field B(t). In principle there are retarda- 
tion effects due to the fact that the dressed molecular state does not change 
instantaneously. Following the above manipulations for time-dependent mag- 
netic field we see that these effects can be neglected if 



h 



d In Z(t) 
dt 



< |e m (t)| 



(149) 



In principle, the Heisenberg equation of motion for the molecular field operator 
also contains an imaginary part due to the fact that the dressed molecule can 
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Fig. 14. Diagrams contributing to the thermodynamic potential of the gas. The non- 
interacting atomic and molecular propagators are denoted by the solid and dashed 
thin lines, respectively. The full molecular propagator is given by the thick dashed 
line. The bare and renormalized atom-molecule coupling constants are denoted by 
the open and filled triangles, respectively. 

decay into a pair of atoms with opposite momenta. The rate for this process 
will be small, however, under the condition given in Eq. (149). We will come 
back to this process when we consider its effect on the coherent atom-molecule 
oscillations. 

For positive detuning the molecular density of states has only a contribution 
for positive energy. For large detuning it is in first approximation given by 

rt ^_ hT m (B)/2 

PmlR '^ vr [(ft, - e k /2 - 6(B)) 2 + (fiTm(S)/2)'] ' 1 J 



where the lifetime of the molecular state is denned by 
m 3. 



r m (B) = 9 ^^{m ■ (i5i) 



As expected, the density of states is, in the case of positive detuning, approxi- 
mately a Lorentzian centered around the detuning with a width related to the 
lifetime of the molecule. It is shown in Fig. 13 by the dashed line. 



4-2 Equilibrium properties 



The equilibrium properties of the gas are determined by the equation of state, 
which relates the total density of the gas to the chemical potential. This equa- 
tion can be calculated in two ways, either by calculating the thermodynamic 
potential and differentiating with respect to the chemical potential, or by di- 
rectly calculating the expectation value of the operator for the total density. 
We discuss both methods, which should, of course, yield the same result. Nev- 
ertheless, to show the equivalence is a subtle matter. 



48 



First, we calculate thermodynamic potential [102]. Within our approximations 
it is in first instance given by the expression 



fi(/i,T) = -Tr 



ln ( G 0,a) 



ln (G- 1 



(152) 



Here, we recall that Go, a (k, icu n ) is the noninteracting atomic propagator of 
the atoms in Eq. (97). The full molecular propagator is given by 

Gm(k,?U;n) = -itw n + e k /2 + 5(B) ~2fi + HX™(ihw n - e k /2 + 2/i) ' (153) 



with the molecular self-energy given in Eq. (111). The so-called ring diagrams 
that contribute to the thermodynamic potential in our approximation are 
given in Fig. 14. The full molecular propagator is denoted by the thick dashed 
line and the noninteracting molecular propagator is denoted by the thin dashed 
line. The noninteracting atomic propagators are indicated by the thin solid 
lines. The total atomic density is calculated by using the thermodynamic iden- 
tity iV = — <9fi(/x, T)/<9/x, which results in 



n = 



d 



tyV Y n [iou n -(e k -ii)/h 



+^f(^„-e k /2 + 2/i))] . 



2/z 



(154) 



After performing the summation over the Matsubara frequencies in this ex- 
pression, the first term corresponds to the density of an ideal gas of bosons. 
The second term in Eq. (154) is more complicated and should, in principle, 
be dealt with numerically. For negative detuning we can gain physical insight, 
however, by expanding the propagator around its pole at the molecular bind- 
ing energy e m (B). This leads to the approximation 



— In 

dfx 



(5 (-ihw n + e k /2 + 5(B) - 2/i + ft£^ B (i/^ n - e k /2 + 2/i 



-2 



l-(^f)'(^ n - ek /2 + 2/i) 



-ihu n + e k /2 + 5(B) - 2/i + KE^(ihuj n - e k /2 + 2/i) 
2 



(155) 



-ifkJ n + e k /2 + e m (B) - 2/i ' 
where we used the expression for the residue of the pole in Eq. (131). With 
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this approximation the sum over the Matsubara frequencies in Eq. (154) is 
performed easily and leads to the result 



1 



n ■ 



EE 



h[3V k „ 
4 



iuj n - (e k - n)/h 



icu n - (e k /2 + e m (B) - 2{j)/h 



y E i N (^ -») + 2iV(e k /2 + e m (B) - 2/x)] 



(156) 



This important result shows that in equilibrium in the normal state and for 
negative detuning the gas in first approximation behaves as an ideal-gas mix- 
ture of atoms and dressed molecules. The same result is found if we neglect in 
the Heisenberg equations of motion for the atomic and dressed molecular field 
operators in Eq. (148) the interaction terms and calculate the total density in 
equilibrium. 

Instead of calculating the thermodynamic potential and differentiating with 
respect to the chemical potential we can also calculate the total density directly 
by using 



n = -G a (x, r; x, r + ) - 2G m (x, r, x, r + ) 



(157) 



An important difference between directly calculating the density in this man- 
ner and calculating it indirectly from the thermodynamic potential is that we 
should use in Eq. (157) not the noninteracting atomic propagator. Instead, we 
should use an approximation to the atomic propagator that contains the same 
self-energy diagrams as the diagrams shown in Fig. 14. Conversely, in calcu- 
lating the thermodynamic potential with the use of Eq. (152) we should not 
use the full atomic propagator. The reason for this is that if we calculate ring 
diagrams with this propagator we find diagrams which are already contained 
in the ring diagram of the full molecular propagator. The following explicit 
example clarifies this further. 

If we use for the atomic propagator the approximation given diagrammati- 
cally in Fig. 15 (a), the ring diagram that contributes to the thermodynamic 
potential is given in Fig. 15 (b). On the other hand, if we use for the molecular 
propagator the approximation given in Fig. 15 (c) the resulting ring diagram, 
given in Fig. 15 (d), is exactly the same as Fig. 15 (b). Clearly, to avoid dou- 
ble counting problems in the calculation of the thermodynamic potential we 
should take only one of these diagrams into account. However, if we calculate 
the density directly from the atomic and molecular propagators we should use 
both the diagrams given in Fig. 15 (a) and (c). 
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Fig. 15. Examples of approximations for (a) the atomic propagator and (c) the 
molecular propagator. The corresponding ring diagrams that contribute to the ther- 
modynamic potential are shown in (b) and (d), respectively. 

We now argue that by directly calculating the density, again for negative 
detuning, we indeed recover the result in Eq. (156). We first calculate the 
contribution arising from the molecular propagator. It is found to be equal to 



n m = -G m (x, r, x, r + ) = -^tU H G m (k, iu n ) 



If 1 



1 



V k J H(3 - iu n - (ftw -2jj)/h 

= I ^3 / d(tuv)p m (k, u)N{hw - 2fi) . (158) 

Taking into account only the pole in the density of states leads to the result 



n r 



Z{B) J J^Y N ^ /2 + UB) ~ 2/U) • (159) 



At first sight this result seems a factor Z(B) to small to agree with the result 
in Eq. (156). However, we have, in fact, already seen in Eq. (144) that the 
contributions from the atoms to the density results in a term proportional to 
2 — 2Z(B). Taking this into account, the result from the direct calculation 
agrees with the result in Eq. (156) obtained previously. 

A different way for obtaining the factor 2 — 2Z(B) in the atomic density is 
to include the self-energy diagram shown in Fig. 16 in the atomic propagator. 
The corresponding mathematical expression is in first instance given by 

ftE a (k, iw n ) = — — G m(k + q, iuj n+m )G a (k, iu n ) . (160) 

q,n 
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Fig. 16. Self-energy of the atoms. The solid and dashed thick lines correspond to the 
full atomic and molecular propagators, respectively. The filled triangles correspond 
to the renormalized atom-molecule coupling constant. 

To understand the physics of this expression, we note that if we neglect the 
energy and momentum dependence of the molecular propagator we have that 
C m (k, iu n ) ~ —H/5(B). Within this approximation the self-energy is given by 
87rn a a rcs (-B)/i 2 /m, which corresponds precisely to the Feshbach-resonant part 
of the self-consistent Hartree-Fock self-energy of the atoms, as expected from 
the diagram in Fig. 16. 

The full calculation of the expression for the self-energy in Eq. (160) is com- 
plicated due to the fact that we have to use the full atomic and molecular 
propagators, which makes the calculation self-consistent. To illustrate in per- 
turbation theory that we are able to reproduce the result in Eq. (144) let us 
simply take the noninteracting atomic and molecular propagators. The self- 
energy is then given by 



V q iftwn - (e k+q /2 - e q + 5(B) - /i) 

To compare with the two-atom calculation for negative detuning performed in 
the previous section, we must take only one other atom present with momen- 
tum —Hk, and no molecules. The self-energy of the atom with momentum hk 
is then given by 

Ag 2 1 

frZ a (k,iuj n ) = — — — — . (162) 

V ihw n - {5(B) - e k - ji) 



With this self-energy the retarded propagator of the atoms is given by 

G[ +) (k,uj) = T1 — f . (163) 

a 1 ^ hw+-e k -*f[hw+ + e k -5(B)]- 1 1 ' 



It has two poles, one close to ek, and one close to 5(B). The residue of the 
latter is given by 

Z k ~ — 5 , (164) 

V [ 2ek - 5(B)] 2 K ] 
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in agreement with the result in Eq. (143). Moreover, we have that 



]TZ k = 2-2Z(5). 



(165) 



k 



Hence, the total density of the atoms is given by 



(2-2Z(S)) 



£ iV( ek /2 + 5(B) - 2/*) + i £ iV(e k - /i) . (166) 




Together with the molecular density from Eq. (159) that becomes 



71- 



£7V(e k + 5(5)-2/x) , 



(167) 



m — 



V 



k 



the total density 2n m + n a is again equal to the result in Eq. (156) to lowest 
order in the interactions. 



4-3 Applications 

In this section we present results on the properties of the normal state of the 
gas. First, we calculate the density of atoms and molecules as a function of the 
detuning, at a fixed temperature. Second, we calculate the density of atoms 
and molecules, and the temperature of the g function of detuning at 

fixed entropy and total density. This calculation is of interest because it gives 
the outcome of a magnetic-sweep experiment through the Feshbach resonance 
in the adiabatic approximation. Finally, we calculate the critical temperature 
for Bose-Einstein condensation as a function of the detuning, at fixed total 
density. 

4-3.1 Density of atoms and molecules 

As we have seen, the density of the gas is most easily calculated by means 
of Eq. (154). We report all our results as a function of the detuning in units 
of the energy g 4 m 3 /A7r 2 h 6 . The temperature is given in units of the critical 
temperature for Bose-Einstein condensation of an ideal gas of atoms with a 
total density n, i.e., 



T = 



3.31fi 2 n 2 / 3 
mk B 



(168) 
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Fig. 17. Fraction of atoms and fraction of atoms in molecules as a function of 
the detuning at a fixed temperature of 2Tq, for two different total densities. The 
solid line is the exact result that includes all two- atom physics, and particular the 
effects of the nonzero lifetime of the molecule at positive detuning, for a total atomic 
density of n = 10 11 cm~ 3 . The dashed line shows the result for this density if we 
approximate the gas by an ideal-gas mixture of atoms and dressed molecules. The 
dotted line shows the exact result for a total density of n = 10 12 cm -3 . 

We compare the exact results, found from numerically performing the summa- 
tion over Matsubara frequencies in Eq. (154), to the ideal-gas mixture result in 
Eq. (156). However, this last equation was derived for negative detuning where 
there is a real pole in the molecular Green's function. We extend this result to 
positive detuning by describing the molecular gas at positive detuning as an 
ideal gas of molecules with a positive bound-state energy given by the energy 
at which the molecular density of states in Eq. (129) has its maximum. It turns 
out that, for small detuning, this maximum is at h 2 /ma 2 . Furthermore, this 
approximation implies that we ignore the physical effects of the lifetime of the 
molecule, that is nonzero for positive detuning, on the equilibrium properties 
of the gas. These lifetime effects are however included in the exact result in 
Eq. (154). 

In Fig. 17 the results of the calculation of the density as a function of the de- 
tuning are shown, for a temperature of T = 2T . Fig. 17 (a) shows the fraction 
of atoms and Fig. 17 (b) shows the number of atoms in molecules, i.e., twice 
the fraction of molecules, as a function of the detuning. The solid and the 
dotted lines show the exact result for a total atomic density of n = 10 11 cm~ 3 
and n = 10 12 cm -3 , respectively. As expected, for negative detuning most of 
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the atoms in the gas are bound to molecules. At positive detuning, most of the 
atoms are free. Moreover, we find that the width in detuning of this crossover 
regime is approximately equal to the temperature. The solid lines show the 
ideal-gas mixture result for a density of n — 10 11 cm -3 , i.e., the result that 
does not incorporate the effects of the nonzero lifetime of the molecules at 
positive detuning. For negative detuning, we observe that this ideal-gas result 
becomes equal to the exact result. This implies that the pole approximation 
in Eq. (155) is indeed a reasonable approximation sufficiently far from res- 
onance. An important conclusion is therefore that, for sufficiently negative 
detuning, we are allowed to treat the gas as an ideal-gas mixture of atoms 
and dressed molecules with binding energy e m (B). For positive detuning, the 
ideal-gas result differs substantially from the exact result. In particular, for rel- 
atively large detuning, the ideal-gas calculation considerably underestimates 
the number of molecules. The exact result shows that there is, even at rela- 
tively large detuning, a significant fraction of molecules in the gas. This is the 
result of the finite lifetime of the molecules in this case. Physically, this comes 
about because the molecular density of states for positive detuning has signif- 
icant spectral weight at low energies. In equilibrium, this leads to a significant 
fraction of molecules. For even larger positive detuning, the ideal-gas result 
reduces again to the exact result. 



4-3.2 Adiabatic sweep through the resonance 

We now calculate the number of atoms and molecules in the gas during an 
adiabatic sweep in the magnetic field, such that the detuning changes from 
positive to negative. The condition for adiabaticity is that the entropy of the 
gas is constant. The entropy is given by 




(169) 



The total number of atoms is, of course, also constant througout the sweep. As 
we have seen, for sufficiently large absolute values of the detuning, the gas is 
well-described by an ideal-gas approximation. For simplicity, we will therefore 
treat the gas here as an ideal-gas mixture since we are mostly interested in 
the final density of atoms and molecules and the final temperature of the gas 
after the sweep, for which an ideal-gas treatment is sufficient [103]. 

In Fig. 18 the results of the calculation of the fraction of atoms and twice the 
fraction of molecules is presented. The total atomic density is taken equal to 
n = 10 13 cm -3 . The solid lines show the result for an initial temperature of 
T = 2To, and the dashed lines show the result for an initial temperature of 
T = 4T . As we go from positive to negative detuning, most of the atoms in the 
gas are converted to molecules. The range of detuning where the conversion 
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Fig. 18. Fraction of atoms and twice the fraction of molecules as a function of the 
detuning for an adiabatic sweep through the resonance. The total atomic density is 
equal to n = 10 13 cm -3 . The solid lines show the result for an initial temperature 
of T = 2T . The dashed lines show the result for T = 4T . 
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Fig. 19. Temperature of the g function of the detuning for a sweep through 

the resonance from positive to negative detuning, for two inital temperatures. The 
total atomic density is equal to n = 10 13 cm~ 3 . 

takes place is proportional to the initial temperature of the gas, as expected. 

In Fig. 19 the temperature is plotted as a function of the detuning for the 
two initial temperatures T = 2T and T = 4T . The total density is again 
equal to n — 10 13 cm -3 . Clearly, the gas is heated as the detuning is changed 
from positive to negative. This is easily understood, since molecules form as 
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Fig. 20. Critical temperature for Bose-Einstein condensation as a function of de- 
tuning. The total density is equal to n = 10 13 cm~ 3 . The solid line shows the result 
of the exact calculations. The dashed line shows the result of treating the gas as an 
ideal-gas mixture. 



the detuning is changed from positive to negative values, and their binding 
energy is released as kinetic energy into the gas. 



4-3.3 Critical temperature 



Finally, we calculate the critical temperature for Bose-Einstein condensation 
of the atom-molecule mixture, at a fixed total atomic density. The results are 
presented in Fig. 20, for a total density of n — 10 13 cm -3 . The solid line shows 
the exact calculation and the dashed line shows the ideal-gas mixture result. 
For positive detuning and far from resonance, we are essentially dealing with 
an atomic gas. Hence we have in this regime that T B ec = T . For sufficiently 
negative detuning we are dealing with a gas of molecules with twice the atomic 
mass, and hence we have that Tbec — 2~ 5//3 T . The feature in the critical tem- 
perature at zero detuning turns out to be a signature of a true thermodynamic 
phase transition, between a phase with a single Bose-Einstein condensate of 
molecules and a phase containing two Bose-Einstein condensates, one of atoms 
and one of molecules, as was first pointed out by Sachdev [104]. This should 
be contrasted with the situation of an atomic Fermi gas near a Feshbach reso- 
nance, where only a BCS-BEC crossover exists [61]. The calculation of the full 
detuning-temperature phase diagram is work in progress and will be reported 
in a future publication [105]. 
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4-4 Many-body effects on the bound-state energy 



In this section we determine the effects of the atomic gas on the molecular 
binding energy. The first step in an examination of these many-body effects is 
the calculation of the molecular self-energy given in Eq. (104). For simplicity, 
we neglect the energy dependence of the atom-molecule coupling constant 
and the many-body effects on this coupling constant. After subtraction of 
the energy-independent shift, the retarded molecular self-energy that includes 
many-body effects is given by the expression 



fisL +) (K.«) = 



n 9 f dk 



1 + iV(e K / 2 +k-AO + N(e K / 2 -k-A0 



frw+ - e K /2 - 2(e k - fj,') 



2e k 



(170) 



Here, we have treated the atoms in the Hartree-Fock approximation which 
effectively implies that the chemical potential is shifted according to 



H =n- 



8na(B)h 2 n a 



/i - 2T 2B n a , 



m 



(171) 



where n a is the density of the atoms. In this expression for the Hartree-Fock 
self-energy correction to the chemical potential we have neglected the energy- 
dependence of the interactions, which is justified as long as the scattering 
length is much smaller than the thermal deBroglie wavelength of the atoms. 

From now on we restrict ourselves to the regime just above the critical temper- 
ature, where we are able to calculate various properties analytically. Since the 
chemical potential approaches zero from below in this regime, we are allowed 
to approximate the Bose distribution function of the atoms by 

N(x) ~ — . (172) 



Within this approximation, the self-energy of the molecules is given by 



^L +) (K, 



2 r dk 1 1 

9 J (2tt)3 hu+ - e K /2 - 2(e k - /i') /?(e K /4 + 6 k - /*') ' ( } 

and we are allowed to also neglect the square-root term that results from the 
first and last terms in the integrand in Eq. (170), and is due to two- atom 
physics. This integral is performed analytically. For h~u < ex/2 — 2/x' the 
self-energy is real and given by 



58 



e K /2-2n' 
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Fig. 21. Graphical solution of the equation for the molecular bound-state energy. 
The solid line indicates the real part of the molecular self-energy as a function of 
fojj. The dashed and dotted lines indicates the function htu — ex/2 — 6(B) + 2/x for 
different values of the detuning 6(B). For hu < £k/2 — 2//, the value of hio at the 
intersections of the dashed and dotted lines with the solid line corresponds to the 
bound-state energy. For hw > €k/2 — 2// it corresponds to the energy of resonant 
states. 



2 m 3/2 



e K /2 - 2/t' - fiw - yJe K /2 - 2/7 



(174) 



For /leu > e K /2 — 2/x' the self-energy contains an imaginary part and is given 
by 



^L +) (k, 



to 



2 m 3/2 



2g 2 m 



e K /2 - 2/i' + i^^; - e K /2 + 2/x' 



(175) 



To find the energy of the molecular state we have to solve for fuv in the 
equation 



fvjj - e K /2 - 8(B) + 2fi - (K, u) = . 



(176) 



A great deal of insight is gained by the graphic representation of this equa- 
tion which is shown in Fig. 21. The solid line represents the real part of the 
molecular self-energy as a function of the energy fko. The straight dashed and 
dotted lines correspond to hio — e K /2 — S(B) + 2/i, for two different values of 
5(B). From this figure it is clear that there is a real solution, i.e., a true bound 



59 



Fig. 22. Molecular bound-state energy as a function of detuning. The dashed line 
show de molecular bound-state energy in vacuum as a function of detuning. The 
solid line shows the many-body effects on the bound-state energy. 

state, if the detuning is such that 



5(B) < 4T 2i X 



2™3/2 



2g m 



4T 2B n a + e c (K)=5 max (K) .(177) 



irh 3 (3Je K /2 - 2// 



Note that this also implies that the position of the resonance in the magnetic 
field is shifted according to 



1 / 2B 2g 2 m 3 / 2 \ 



(178) 



due to many-body effects. 

For a magnetic field such that the detuning is just below the maximum value 
^max(K) given in Eq. (177), the bound-state energy is given by 



frw 



K 



-2// 



V CO 



+ 



h 2 K 2 

2m cS 



(179) 



with an effective mass given by 
~3(5(B) - 4T 2B n a 



m cS = 2m 



e r (0) 



1 - 



2 (5(B) - AT n 

3 TJo) 



(180) 



This effective mass has a minimum value of 4m/3 at detuning 5(B) = 4T 2B n a + 
3e c (0)/4, and diverges for smaller detunings close to 4T 2B n a . In the limit of the 



60 



Pm ( k > ro ) 




Z 8 (fito-Fito K ) 



A 



M0 K 



FiCO + 



Fig. 23. Molecular density of states with many-body effects. Apart from the delta 
function that corresponds to the bound state there are two resonant states, indicated 
by the dashed lines. 

detuning 5(B) — > — oo we have to recover the two-body bound state with mass 
2m, which shows that this divergence is due to the approximations we have 
adopted. As already discussed, we have in particular neglected the first and 
last terms in the integrand in Eq. (170) that result from two-atom physics. 
Nevertheless, the fact that the effective mass is smaller than the mass of 
a molecule close to resonance indicates that the molecule crosses over to a 
more complex many-body bound state of the system. Precisely at the shifted 
resonance at 5(B) = AT 2B n a + e c (0) the effective mass is again equal to 2m. 
Another interesting feature of the excitation is that for a given detuning it 
only exist at small momenta such that Eq. (177) is obeyed. 

The intersections at energies hw > e K /2 — 2/jf in Fig. 21, as for example shown 
by the dotted line, correspond to resonant states since the self-energy contains 
an imaginary part at these energies. The energies of these resonant states is 
determined by solving for hw in the equation 



hw - e K /2 + 2/JL- 6(B) + 



2f7 2 m 3/2 ^ eK /2 - 2/i' 

Trh 3 P hw 



= . 



(181) 



For a detuning that obeys the condition in Eq. (177) and such that 



5(B) > 2/i - e K /2 + 



8g 2 m 3 / 2 



e K /2 - 2/x' = <5 min (K) , 



(182) 



\ irh 3 (3 



there are two solutions of this equation. They are given by 
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tuv± = -(e K /2 + 6(B)-2 f i) 

8g 2 m 3/2 - 2yl \ 

tt&P (e K /2 + 6(B)-2^ ' 1 J 

/ 

For large detuning we have that fru + ~ e K /2 + 6(B) — 2/i, from which we see 
that this resonant state physically corresponds to the bare molecular state, 
which has obtained a finite lifetime due to the interaction with the atomic 
continuum. The resonant state at energy ftw_ is not present in the two-atom 
case but arises purely due to many-body effects. This situation is somewhat 
similar to the Kondo-resonant state that arises in a Fermi gas near a Feshbach 
resonance [68]. 

An illustration of the many-body effects on the molecular bound-state en- 
ergy is shown in Fig. 22. The dashed line indicates the situation in vacuum. 
For negative detuning there is a true molecular state whose energy depends 
quadratically on the detuning, as given in Eq. (126). For positive detuning the 
molecule has a finite lifetime and therefore corresponds to a resonant state, 
whose energy is in first approximation equal to the detuning. Due to many- 
body effects, the position of the Feshbach resonance is shifted. Nevertheless, 
there is still a molecular state with an energy dependence that is quadratic on 
the many-body renormalized detuning. However, for a detuning larger than 
6 m i n but less than <5 max this molecular state coexists with two resonant states, 
one close to the detuning and one just above the continuum threshold. The 
molecular density of states for the latter situation is shown in Fig. 23. The 
delta function corresponds to the molecular bound state. The dashed lines 
indicate the position of the resonances. For large positive and large negative 
detuning the many-body effects are negligible and the result reduces to the 
two- atom result. 

Finally, we remark that the resonant state at energy hjj_ , that arises solely due 
to many-body effects, leads to a nonzero number of bare molecules, even if the 
temperature is much smaller than the detuning. This effect can be measured by 
directly measuring the number of bare molecules, as achieved recently by Chin 
et al. [106]. The investigation of the magnitude and temperature dependence 
of this effect is intended for future work. 



5 Mean-field theories for the Bose-Einstein condensed phase 

In the first part of this section we derive the mean-field theory that results 
from our effective quantum field theory. This mean-field theory is appropriate 
for the description of the Bose-Einstein condensed phase of the gas. In the 
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second section we discuss other possible mean-field theories and discuss the 
similarities and differences between them and our mean-field theory. 



5.1 Popov theory 



In this section we derive the mean-field equations for the atomic and molecular 
condensate wave functions. In the first part of this section we derive the time- 
independent equations and discuss the excitation spectrum. In the second part 
we derive the time-dependent mean-field equations. 



5.1.1 Time-independent mean-field equations 

The mean-field equations for the atomic and molecular condensate wave func- 
tions are derived most easily by varying the effective action in Eq. (112) with 
respect to a£ n and &£ n , respectively. Before doing so, however, we remark 
that an important property of this effective action is its invariance under 
global U(l) transformations. Namely, any transformation of the form 



b Kn ^b Kn e 2id , (184) 

with 8 a real parameter, leaves the action unchanged. The conserved quan- 
tity, the so-called Noether charge, associated with this invariance is the total 
number of atoms. The appearance of the atomic and the molecular conden- 
sates breaks the £7(1) invariance since the wave functions of these condensates 
have a certain phase. According to Goldstone's theorem, an exact property of 
a system with a broken continuous symmetry is that its excitation spectrum 
is gapless [107]. Since our mean-field theory is derived by varying a U{1)- 
invariant action, this property is automically incorporated in the mean-field 
theory. 

To derive the time-independent mean-field equations, that describe the equilib- 
rium values of the atomic and molecular condensate wave functions, we substi- 
tute into the effective action ao,o — > (fi&V P^V + a o,o an d &o,o — ¥ 4>m \/]3ftV + &o,o- 
Here, a and <f) m correspond to the atomic and molecular condensate wave 
functions, respectively. Requiring that the terms linear in a ,o an d 6 ,o vanish 
from the effective action leads to the equations 



/i0 a =T b 2 g B (2/i - 2ft£ HF ) |0 a | 2 0a + 2 [g 2B (2/i - 2ft£ HF )] * 0> m , 
2/i0 m =[<5(5) + hY™ (2/i - 2/ffi HF )] m + g 2B (2/i - 2/lS HF ) 2 . (185) 
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Fig. 24. Hartree-Fock self energy of the atoms. The dotted lines correspond to 
condensate atoms. The dashed line corresponds to the full molecular propagator. 

A crucial ingredient in these equations is the Hartree-Fock self-energy of the 
noncondensed atoms. This self-energy is the mean-field energy felt by the 
noncondensed atoms due to the presence of the atomic condensate. Taking 
into account the energy- dependence of the interactions, it is determined by 
the expression 

with n a = 1 a | 2 the density of the atomic condensate. Its diagrammatic rep- 
resentation is given in Fig. 24. The overall factor of two comes from the con- 
structive interference of the direct and exchange contributions. Far off reso- 
nance we are allowed to neglect the energy-dependence of the effective atom- 
atom interactions, and the Hartree-Fock self-energy of the atoms is given by 
8na(B)ti 2 n Si /m, as expected. The Hartree-Fock self-energy is essential for a 
correct description of the equilibrium properties of the system. The physical 
reason for this is understood as follows. In the condensed phase the chemi- 
cal potential is positive. The energy of a condensate molecule is equal to 2fi, 
which is therefore larger than the continuum threshold of two atoms in vac- 
uum. Without the incorporation of the Hartree-Fock self-energy, the molecular 
condensate would therefore always decay and an equilibrium solution of the 
mean-field equations would not exist. However, due to the presence of the 
atomic condensate the continuum threshold shifts by an amount 2/i£ HF , and 
the molecular condensate is stable. 

To study the collective excitation spectrum over the ground state determined 
by Eq. (185), we consider the effective action up to second order in the fluc- 
tuations, which is known as the Bogliubov approximation [108]. To facilitate 
the notation we introduce the vector u k n by means of 
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Uk,n = 



( n \ 



(187) 



With this definition, the quadratic part of the effective action is given by 



5 B [u f ,u] = --J2 u tn ■ G B \k,iu n ) u k , n , 



k,n 



(188) 



where the Green's function of the fluctuations is determined by 



c _1 c _1 

a coup 
^coup u m 



(189) 



The atomic part of this Green's function is found from 



-^G a 1 (k, iu;„ 



^G 0)a (k, iu n ) 



\ 







+ 



2T 2B (i/^ n 

->2B 



e k /2 + 2/i') n a T b 2 g B (2/.') </> 2 + 2 ^ 2B (2/x' 



^ 2B (2//) (0:) 2 + 2^ 2B (2/i')C 2T 2B (ihw n - 6 k /2 + 2/x') n 



(190) 



where // = // — /i£ HF . Note that in the absence of the coupling to the molec- 
ular condensate, this result reduces to the well-known result for the Green's 
function that describes phonon propagation in a weakly-interacting Bose con- 
densate. We have in this case, however, also explicitly taken into account 
the energy dependence of the coupling constants. Therefore we know that in 
the limit of vanishing coupling g 2B the propagator in Eq. (190) has a pole 
that determines the gapless dispersion relation for the phonons. For energy- 
independent interactions this so-called Bogoliubov dispersion is given by 



/ 87ra bg ^ 2 n a 

'Ck "1 6 k • 



m 



(191) 



The molecular part of the Green's function Ge(k, iu n ) is determined by 
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G m 1 (k,io; n ) 



\ 



Gm^k, -iw n ) ) 
where the single-molecule propagator is given by 



(192) 



- hG m \k, iu n ) = -ifrw n + e k /2 + 5(B) - 2fi 



+KT? B (ihw n - e k /2 + 2/i - 2ft£ HF ) 



(193) 



From the previous section we know that the Green's function in Eq. (193) for 
negative detuning has a pole at the molecular binding energy. There are now, 
however, mean-field effects on this binding energy due to the presence of the 
atomic condensate, incorporated by the Hartree-Fock self-energy ft£ HF [85]. 
Finally, the Green's function that describes the coupling between the atomic 
and molecular fluctuations is given by 



-tiG*(k,iu n ) 



g 2B (ihLu n -e k /2 + 2fi')\ ^ 

2g 2B (i?ku n -e k /2 + 2n , ) ( j) a J 



(194) 



The spectrum of the collective excitations is determined by the poles of the 
retarded Green's function for the fluctuations Ge(k, u + ). This implies that 
we have to solve for hu in the equation 

detGi^k,^) = . (195) 



This is achieved numerically in the next section to determine the frequency of 
the Josephson oscillations between the atomic and the molecular condensate. 
However, we are already able to infer some general features of the excitation 
spectrum of the collective modes. We have seen that in the absence of the 
coupling between the atomic and molecular condensate, we have that one dis- 
persion is equal to the gapless Bogoliubov dispersion with scattering length 
dbg. In the presence of the coupling this branch corresponds again to phonons, 
but the dispersion is now approximately equal to the Bogoliubov dispersion 
for the full scattering length a(B). There is a second dispersion branch that 
for small coupling g 2B lies close to the molecular binding energy. At nonzero 
coupling this branch corresponds to coherent atom-molecule oscillations, i.e., 
pairs of atoms oscillating back and forth between the atomic and molecular 
condensate. Physically, the difference between the two branches is understood 
by realizing that for the phonon modes the phases of the atomic and the 
molecular condensate are locked to each other and oscillate in phase. Since the 
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action is invariant under the transformations in Eq. (184) we conclude that the 
phonons are indeed gapless, and, in fact, correspond to the Goldstone mode 
associated with the breaking of the U(l) symmetry by the condensates. For 
the coherent atom-molecule oscillations the phases of the atomic and molec- 
ular condensate oscillate out of phase and hence the associated dispersion is 
gapped. As a final remark we note that we indeed have that 

detGi^O) = , (196) 



which shows that there is indeed a gapless excitation, in agreement with Gold- 
stone's theorem. 



5.1.2 Time- dependent mean- field equations 

The time-dependent mean-field equations are found most easily by taking the 
expectation value of the Heisenberg equations of motion in Eq. (117). For 
notational convenience we restrict ourselves to the situation that we are close 
to resonance and hence neglect the energy-dependence of the various couplings. 
Moreover, we only take into account the leading-order energy dependence of 
the molecular self-energy, as given in Eq. (116). Furthermore, we assume that 
we are at such low temperatures that the effects of the thermal cloud may be 
neglected. Within these approximations, the mean-field equations are given 
by 



<90 a (x,t) 



h 2 V 2 4na hK h 2 . , . . ., 

+ M0 a (x,t)| 2 



2m m 
+ 2 #I(x,t)0 m (x,t) , 
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ih- 



(M) 
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h 2 v 2 

4m 



+ 6(B(t)) 



&n(M) + 00a(M) 
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2 m 3 / 2 Id h 2 v 2 
%\ih— + 



2nh 3 



dt 



Am 



2/l£ HF m (x,t) . 



(197) 



Note that, since we use renormalized coupling constants in these equations, 
we should not explicitly include also the so-called anomalous averages because 
this leads to double-counting of the interatomic interactions. This is explained 
in detail in the next section. 

The equilibrium solutions of these mean-field equations are space-independent 
and of the form 



a (x,t) = a e 

> m (x,t) = m e" 



-ifit/h 
-2i/it/h 



(198) 
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Substitution in Eq. (197) reproduces the time-independent equations for a 
and m within the above approximations. Moreover, by linearizing around 
these equilibrium solutions we find again the collective-mode spectrum dis- 
cussed in the previous subsection. 

We now discuss the solution of the homogeneous version of the time-dependent 
mean-field equations in Eq. (197). These equations are given by 



ih- 



dt 



S(B(t))-g 



2nh 3 



dt 



2ft£ HF 



d4> a (t) _4ira hg h 2 



ih 



dt 



m 



W| 2 0aW+2#:(t)0 m (t) 



(199) 



Two different situations can occur, that of time-independent detuning and that 
of time-dependent detuning. Let us first discuss the case of time-independent 
detuning. In this case we are able to solve the equation for the molecular 
condensate wave function by introducing the Fourier transform of the zero- 
momentum part of the retarded molecular Green's function. This Fourier 
transform is, for the most interesting case of negative detuning, given by 



Gf)(t-f) = /^GW(0 )W )e^) 

i9(t-t')g 2 m*/ 2 °fdu v / ^e-< w+2EHF ) (t - t ' ) 



7T 



h 2 



2% [hoj + 2ft£ HF - 5(B)} 2 + (# 4 m 3 /47r 2 ft 6 )^ 



-i9(t-t')Z(B)exp --e m (B)(t-t') 

h 



(200) 



where e m (B) is the molecular binding energy that includes also the effects of 
the Hartree-Fock self-energy. The molecular condensate wave function is, in 
terms of this Green's function, given by 



Mt) = f / dt'G^\t - f )#(f ) + m (O)e— ^ , (201) 
o 

for t > 0. This result is substituted in the equation for the atomic condensate 
wave function, which can subsequently be solved numerically. 

The second situation we can have is that of a time-dependent detuning. To take 
into account the fractional derivative acting on the molecular wave function 
in the second equation in Eq. (199), we use its definition in frequency space. 
Hence we have that 
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— oo — oo 
oo oo 

= I dt' I ^Vf^e-^'^Ut') ■ (202) 

— oo — oo 

This specific definition is referred to in the literature as the Weyl definition of 
a fractional derivative [109]. Unfortunately, the integral over uj in the above 
expression does not converge. This problem is overcome by considering also the 
next-order energy-dependence of the molecular self-energy. Therefore, we take 
for the molecular self-energy the expression in Eq. (124), i.e., the molecular 
self-energy with the effective range r\ g = 0. The equation for the molecular 
mean field is then given by 



d 



"dt 



Hi— - 8(B(t)) + - 
dt l_,K^^!_ 2 /ffiHF 



0m(t) = g<f>l(t) . (203) 



The term that involves the fractional derivatives is now rewritten as 



,/ 0m(t) 

1 _ i K s h^^ ih a__ 2hE HF 

= 7 dt , [ ^ i^^-2f^e-^M^) m) 
_7 J 2tt 1 _,Kglv^ v /^_ 2 ^ S HF 

For large u the integrand becomes equal to a constant which gives rise to a 
delta function 5(t — t'). Taking this into account, the final result for this term 
is given by 
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(205) 



where the characteristic time r = ma1 g /h and the complementary error func- 
tion is defined by means of 



oo 

Erfc(z) = -= / dwe"™ 2 = 1 - Erf (2) . 



(206) 
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This final result shows that the term involving the fractional derivatives may 
be dealt with numerically as a term that is nonlocal in time. In the next 
section we present results of numerical solutions of the time-dependent mean- 
field equations using the Green's function method. 



5.2 Hartree-Fock-Bogoliubov theory 



A completely different approach to arrive at mean-field equations that describe 
the Bose-Einstein condensed phase of a system with Feshbach-resonant inter- 
actions has been put forward by Kokkelmans and Holland [79] and Mackie et 
al. [78]. Their treatments are physically similar but differ in some technical 
details. We discuss here the approach of Kokkelmans and Holland. 

Their starting point is the microscopic atom-molecule hamiltonian in Eq. (87). 
The first step is to approximate the interatomic potential and the atom- 
molecule coupling as contact interactions, according to 



\/ TT (x-x')~W(x-x') , 

g u (x-x!)~go8(x-x!) . (207) 

Roughly speaking, this approximation is validated by the fact that the de- 
Broglie wavelength of the atoms and molecules is much larger than the range 
of the interactions. However, the use of contact interactions leads to ultraviolet 
divergencies in the theory which have to be regularized by introducing a ultra- 
violet cut-off k\ in momentum space. The unknown microscopic interaction 
parameters Vq and go are then expressed in terms of the experimentally known 
parameters g, A/x, and a^ g , and the cut-off in such a way that the final 
equations correctly describe the two-atom physics and are cut-off independent 
in the limit of a large cut-off. This renormalization procedure is discussed in 
detail below. 



First we derive the so-called Hartree-Fock-Bogoliubov equations of motion. 
Within the above approximation, the hamiltonian for the system is given by 
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h 2 v 2 v ?y 
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rfx vi(x 

+ J dx?/v 
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^a(x) 



4m 



+ 1/(5) 



^m(x) 



(208) 



where u(B) is a bare and also cut-off dependent detuning for the molecular 
state. In this hamiltonian, the Schrodinger operators that annihilate an atom 
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and a molecule are denoted by V> a (x) and ^ m (x), respectively. Their hermitian 
conjugates are the creation operators. 

The starting point in the derivation of the Hartree-Fock-Bogoliubov equations 
of motion are the equations of motion for the Heisenberg operators ^> a (x, t) 
and ^ m (x, t), that follow from the hamiltonian in Eq. (208). They are given 
by 



%h dt 

..dxj> m (x,t) 



h 2 v 2 

' 2m 

h 2 v 2 

4m 



+ Wl(x,t)4CM) 



^ a (x,t) +2^1(x,t)^ m (x,t) , 



+ u(B) 



(209) 



The next step is to separate out the expectation value of the Heisenberg op- 
erators. These expectation values are constant in space since we are dealing 
with a homogeneous system. We write the Heisenberg operators as a sum of 
their expectation values and an operator for the fluctuations according to 



V>a(x, t) = (V> a (x, t)) + Xa(x, t) = (f) a (t) + Xa(x, t) , 

^ m (x, t) = (^ m (x, t)) + Xm(x, f) = m (t) + Xm(x, *) . (210) 

We substitute this result into the Heisenberg equations of motion and take 
the expectation values of these equations. These expectation values are then 
decoupled in a manner that is similar to Wick theorem. This is, of course, an 
approximation in this case since we are dealing with an interacting system. In 
detail, we only take into account the expectation values {ip a ), (0m), (%aXa), 
and (xlxa)- This leads to four coupled equations of motion for these expecta- 
tion values. We define the so-called normal and anomalous expectation values 
according to 



G N (r,t) = (xi(x,t)x a (x , ,t)) , 

G A (r,t) = (Xa(x,t)x a (x , ,t)) , (211) 

which only depend on the difference r = x — x' due to translational invari- 
ance of the system. Note that the normal average yields the density of non- 
condensed atoms according to n'(t) = G^(0,t). Including the normal average 
does not alter the conclusions of the following discussion. Therefore, we assume 
from now on that we are at such low temperatures that there is essentially no 
thermal cloud present, and therefore take G^(r,t) = 0. 

The Hartree-Fock-Bogoliubov equations of motion are given by 
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(212) 



Note that, as they stand, these equations cannot be derived by varying a 
[/(l)-invariant action. However, we have seen that this U(l) invariance is an 
exact property of the theory. This problem is overcome by realizing that the 
anomalous average G A is in fact proportional to the atomic condensate wave 
function, since it is zero in the normal phase of the gas. More precisely, we 
have that G A oc 4> 2 which renders the equations for the atomic and molec- 
ular condensate wave function L^l) -invariant. Moreover, elimination of the 
anomalous average for the Hartree-Fock-Bogoliubov equations of motion in 
Eq. (212) leads to renormalization of the bare couplings V and g . We have 
already seen in Section 3.1 that introducing a pairing field into the theory 
leads to a summation of the ladder Feynman diagrams. We expect something 
similar to occur in this case [110,111]. 

To study how this renormalization works in detail we study the equilibrium 
solutions of the Hartree-Fock-Bogliubov equations. Therefore, we substitute 



m (t) = m e- 2i ^, 
G A (r,t) = G A (r)e- 2i ^ h , 



(213) 



from which we find the time-independent Hartree-Fock-Bogoliubov equations 



W>a = WalVa + [V O G A (0) + 2<?o0 m ] 4>l , 

2/Mf) m = is(B)(f) m + g [4>\ + G A (0)] , 

h 2 v 2 



2/2G A (r) 



m 



4W a 



Ga(t) 



+ 



V o( f ) l + V o G A (O)+2g o m 5(r) . 



(214) 



The equation for the anomalous average G A (r) is solved by Fourier transfor- 
mation. This gives the result 



G A (0) 



Vo 

V 



^l k l< fc A 2 M +-2e k -4V o |0a| 2 



V ^|k|<fc A 2/j+-2e k -4Vb|</>a| 2 



^0 



(215) 
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which explicitly shows that the anomalous average is proportional to the 
atomic condensate wave function. Note also that we have to regularize this 
expression by using the ultraviolet cut-off k^, since it would be ultraviolet 
divergent otherwise. Converting the sum over momenta to an integral, we find 
the final result for the anomalous average 



Gm = ^ + V-' ^ (216) 



Substitution of this result into the equations of motion for the atomic and 
molecular condensate wave functions gives in first instance 

2n4> m = v T (B)4> m + g T 4>l , (217) 
where the renormalized interaction and atom-molecule coupling are given by 



# r = v 3/2 / — + #o , 218 

and the renormalized detuning is given by 



MB) = 2* f + 



Finally, we have to express these renormalized quantities in terms of the ex- 
perimentally known parameters ab g , g, and 5(B). Moreover, this has to be per- 
formed in a manner that does not depend on the cut-off in the limit k^ — > oo. 

The renormalization procedure used by Kokkelmans and Holland is given by 



i _ mk A 4:TTa bg h 2 ' 

9 

^° y — mfc A ina hg h 2 ' 
27T 2 h 2 m 



m = m + ( 220 ) 
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Eliminating the microscopic parameters Vo, go, and v(B) in favor of a bg , g, 
and 5(B) finally yields the renormalized mean- field equations for the atomic 
and molecular wave functions 



W>a = 
2/i0 m = 



g2 m 3/2 



1 + (2/i-4\/ O |0 a | 2 

5(3) 
-1 



^a| 2 0a + 



2<? 



1 + 

v /2/i-4V o |0 a | 2 



f (2/i-4\/ o |0 a | 2 ) 



aVm ; 



2^ l-^|a bg | v /f(2 /U -4V o |0 a | 2 )_ 

5 ,2 



l + m bg ^/f(2/x-4Vo|0 a | 2 ) a ' 

where we have retained the term 4Vo|0 a | 2 in the energy arguments of the 
coupling constants. In the limit k\ — > oo this term vanishes and the above 
renormalized equations no longer depend on the microscopic parameters and 
the cut-off. 



(221) 



The above equations are very similar to the mean-field equations of our ef- 
fective field theory in Eq. (185), if we neglect the effective range of the in- 
teractions in the couplings and the self-energy of the molecules in the latter 
equations. There is, however, another and much more important difference 
between the two mean-field theories. In the mean-field theory that we have 
derived from our effective quantum field theory we have included the Hartree- 
Fock self-energy that is due to the mean-field interactions of the condensate 
on the thermal atoms. This Hartree-Fock self-energy is crucial for a correct 
description of the equilibrium properties of the system. In the Hartree-Fock- 
Bogoliubov equations the Hartree-Fock self-energy is replaced by the energy 
4Vo|0 a | 2 , which corresponds to the mean- field energy resulting from the un- 
renormalized interaction. The fact that the interaction between the condensed 
and noncondensed atoms is not renormalized is a well-known problem of the 
Hartree-Fock-Bogoliubov theory [112]. Note also that for a nonzero effective 
range r bg the two-atom physics is not incorporated exactly, and this will lead to 
a discrepancy with experiment as shown in the following section. Although the 
renormalization of the interactions between condensate atoms is, for r bg = 0, 
correctly achieved, the interactions between condensate atoms and thermal 
atoms is not correctly incorporated. In the limit where the cut-off k\ goes 
to infinity this mean-field energy actually vanishes and we conclude from our 
previous discussion that the Hartree-Fock-Bogoliubov equations in Eq. (214) 
have no equilibrium solution. As a result also a linear-response analysis, sim- 
ilar to the one carried out in Section 6, is not possible with this approach. 
Moreover, the above renormalization procedure relies on the presence of the 
anomalous average Ga(t) which makes the theory inapplicable above the crit- 
ical temperature for Bose-Einstein condensation. Hence also a description of 
the thermal cloud of a Bose-Einstein condensed gas cannot be obtained in this 
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manner. Note also that the above result explicitly shows that the inclusion of 
the pairing field G\(r) indeed leads to the summation of the ladder diagrams. 
This is the reason why it is exact not to include anomalous averages in our 
mean-field equations. Their effect is already incorporated by using properly 
renomalized coupling constants. 

Finally, we make some remarks about the theory put forward by Kohler et 
al. [84]. These authors do not explicitly include the molecular field respon- 
sible for the Feshbach resonance into their theory, but instead use a sep- 
arable pseudopotential for the interaction between the atoms that, when in- 
serted in the Lippmann-Schwinger equation, reproduces the energy- dependent 
T-matrix. Subsequently, they use the single-channel version of the above- 
described Hartree-Fock-Bogoliubov theory to arrive at their mean-field equa- 
tions. The theory of Kohler et al. is derived from our effective atom-molecule 
approach by neglecting the effect of the molecular condensate on the atoms. 
The molecular field can then be integrated out, which leads to an energy- 
dependent T-matrix for the atoms. We have seen in Eq. (128) that close to 
resonance the energy-dependence of this T-matrix is equivalent to the energy- 
dependence of the T-matrix in the single-channel case. Close to resonance, 
therefore, the mean-field theory of Kohler et al. incorporates the correct two- 
atom physics. However, their approach cannot fully recover all the properties 
of the molecules, which have been integrated out of the problem. This can for 
instance be seen from the fact that the theory contains only the ratio g 2 /Afi 
instead of the independent quantities g and A/i, seperately. Their theory also 
does not incorporate the mean-field shift on the noncondensed atoms due to 
the atomic condensate, as we have seen explicitly above. The latter feature 
again disables a linear-response analysis of the beautiful experiments we are 
going to discuss next. 



6 Coherent atom-molecule oscillations 



In this section we discuss the experimental observation of atom-molecule co- 
herence in a Bose-Einstein condensate [69,80], and its theoretical description 
in terms of the mean-field theory derived in the previous section. In the first 
section we discuss the experimental results. In the next section we calculate 
the magnetic-field dependence of the frequency of the coherent atom-molecule 
oscillations in linear-response theory. In the final section we present the results 
of calculations that go beyond this linear approximation. 
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Fig. 25. Typical magnetic-field pulse sequence as used in the experiments of Donley 
et al. [69] and Claussen et al. [80]. 

6. 1 Experiments 



In the experiments of Donley et al. [69] and Claussen et al. [80], performed 
both in Wieman's group at JILA, one makes use of the Feshbach resonance 
at B = 155.041(18) G(auss) in the [/ = 2;m f = -2) hyperfme state of 85 Rb. 
The width of this resonance is equal to AB = 11.0(4) G and the off-resonant 
background scattering length is given by a^ g = — 443ao, with ao the Bohr 
radius. The difference in the magnetic moment between the open channel and 
the closed channel is given by A/x = — 2.23/ib, with /i B the Bohr magneton 
[79]. 

In both experiments, one starts from a stable and essentially pure conden- 
sate of about iV c = 10000 atoms at a magnetic field such that the effective 
scattering length is close to zero. This implies that, since the condensate is 
in the noninteracting limit, its density profile is determined by the harmonic- 
oscillator groundstate wave function. The harmonic external trapping poten- 
tial is axially symmetric, with trapping frequencies v r = 17.4 Hz and v z = 6.8 
Hz in the radial and axial direction, respectively. 

Starting from this situation, one quickly ramps the magnetic field to a value 
-Bhoid close to the resonant value and keeps it there for a short time thoid before 
ramping to a value -B ev oive- The magnetic field is kept at this last value for a 
time t C voivc before performing a similar pulse to go back to the initial situation. 
The duration of all four magnetic- field ramps is given by t ramp . A typical pulse 
is illustrated in Fig. 25. Both the ramp time t ramp and the hold time thoid are 
kept fixed at values of 10—15 /is. The time t e voive between the pulses is variable. 

Such a double-pulse experiment is generally called a Ramsey experiment. Its 
significance is most easily understood from a simple system of two coupled 



76 



harmonic oscillators. Consider therefore the hamiltonian 



H=\(tf St). 




(222) 



where a) and ¥ create a quantum in the oscillators a and b, respectively, and 
A denotes the coupling between the two oscillators. 

We consider first the situation that the detuning 5(t) is time independent. The 
exact solution is found easily by diagonalizing the hamiltonian. We assume 
that initially there are only quanta in oscillator a and none in b, so that we 
have that (6*6) (0) = 0. The number of quanta in oscillator a as a function of 
time is then given by 



(0*0) (t) 



(hw) : 



■sin 2 (wt/2) 



(a t a)(0) , 



(223) 



with the frequency w given by 
hw = VS 2 + A 2 . 



(224) 



We see that the number of quanta in the oscillator a oscillates in time with fre- 
quency zu. Such oscillations are called Rabi oscillations. Note that the number 
of quanta in oscillator b is determined by 

ftb)(t) = --^- 2 sin 2 (wt/2) (ata>(0) , (225) 



so that the total number of quanta is indeed conserved. 

Suppose now that we start from the situation with all quanta in the oscillator 
a and none in b and that the detuning is such that 8(t) ^> A. Then we have 
from Eq. (223) that (tfa)(t) ~ (a t a)(0) and (6*6) (t) ~ 0. Starting from this 
situation, we change the detuning instantaneously to a value 5(t) ~ and 
keep it at this value for a time thoid- During this hold time quanta in oscillator 
a will go to oscillator b. Moreover, if t ho i d is such that 

(226) 



on average half of the quanta in oscillator a will go to oscillator b. Such a pulse 
is called a 7r/2-pulse. The defining property of a 7r/2-pulse is that it creates 
a superposition of the oscillators a and b, such that the probabilities to be in 
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oscillators a and b are equal, and therefore equal to 1/2. This is indicated by 
the average (a*b)(t). At t — this average is equal to zero because there is 
no superposition at that time. We can show that after the above 7r/2-pulse 
the average (atb)(t) reaches its maximum value. In detail, the state after the 
7r/2-pulse is equal to 



1 

W\ 



"at +£t 

~7T 



-, N 



10), 



(227) 



where the ground state is denoted by |0), and N = (a^a)(Q) . 

We can now imagine the following experiment. Starting from the situation 
5(t) ^> A, we perform a 7r/2-pulse. Then jump to a certain value <5 evo i ve for 
a time tcvoivc 5 and after this perform another 7r/2-pulse and jump back to 
the initial situation. The number of quanta in the oscillator a, a measurable 
quantity, then oscillates as a function of t C voivc with the oscillation frequency 
determined by Eq. (224) evaluated at the detuning <5 ev oivc- The second it/ 2- 
pulse enhances the contrast of the measurement thus providing a method of 
measuring the frequency w as a function of the detuning with high precision. 

This is basically the idea of the Ramsey experiments performed by Donley 
et al. [69] and Claussen et al. [80]. Roughly speaking, the atomic conden- 
sate corresponds to oscillator a and the molecular condensate to oscillator b. 
Therefore, after performing the double-pulse sequence in the magnetic field 
one makes a light-absorption image of the atomic density from which one ex- 
tracts the number of condensate and noncondensed atoms. Since this imaging 
technique is sensitive to a specific absorption line of the atoms it does not 
measure the number of molecules. 

From the above discussion we expect to observe oscillations in the number of 
condensate atoms. Moreover, if the situation is such that the detuning between 
the pulses is relatively large the effect of the coupling can be neglected and 
the frequency of the observed oscillations corresponds to the energy difference 
between the atoms and the molecules, i.e., the molecular binding energy. This 
is indeed what is observed, thereby providing compelling evidence for the 
existence of coherence between atoms and molecules. 

In Fig. 26 the experimental results of Claussen et al. [80] are presented. Fig. 26 
(a) and (b) show the number of atoms in the atomic Bose-Einstein conden- 
sate as a function of t cvo i ve after a double-pulse sequence. Clearly, there is an 
oscillation in the number of atoms in both cases. In Fig. 26 (a) the magnetic 
field between the pulses is -B ev oivc = 156.840(25) G. In Fig. 26 (b) we have 
-^evolve = 159.527(19) G which is further from resonance. This explains also 
the increase in frequency from (a) to (b) since further from resonance the 
molecular binding energy is larger. 
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Fig. 26. Experimental observation of coherent atom-molecule oscillations. The fig- 
ures are taken from Ref. [80]. Figures (a) and (b) show the number of atoms in 
the atomic condensate as a function of the time between the two pulses in the 
magnetic field. The solid line indicates the fit in Eq. (228). For (a) we have that 
-^evolve = 156.840(25) G. The frequency and damping rates are respectively given 
by v e = 2tt x 0.58(12) kHz, a = 7.9(4) atom/^s, and = 2tt x 0.58(12) kHz. For (b) 
the magnetic field £ evo i ve = 159.527(19) G and u c = 157.8(17) kHz. The damping 
is negligible for the time that is used to determine the frequency. Note that the 
frequency has increased for the magnetic field further from resonance. Figures (c) 
and (d) show the observed frequency of the coherent atom-molecule oscillations as a 
function of the magnetic field. The solid line is the result for the molecular binding 
energy found from a two-body coupled-channels calculation using the experimental 
results for the frequency to accurately determine the interatomic potential [80]. Only 
the black points were included in the fit. The inset shows that, close to resonance, 
the observed frequency deviates from the two-body result. 

What is also observed is that there is a damping of the oscillations and an 
overall loss of condensate atoms. Experimentally, the number of atoms in the 
condensate is fit to the formula 

N c (t) = iV avcragc - at + A exp(-/ft) sin(cu c t + 0) , (228) 

where iV avera ge is the average number of condensate atoms, A and <fi are the 
oscillation amplitude and phase, respectively, and (5 is the damping rate of 
the oscillations. The overall atom loss is characterized by a rate constant a. 
The experimentally observed frequency is equal to u E = lnJvl — [l3/2n] 2 . By 
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defining the frequency of the coherent atom-molecule oscillation in this way 
one compensates for the effects of the damping on the frequency. For the results 
in Fig. 26 (a) we have that (3 = 2tc x 0.58(12) kHz and a = 7.9(4) atom///s. 
The frequency is equal to u c = 9.77(12) kHz. For Fig. 26 (b) the frequency is 
equal to u e = 157.8(17) kHz. The damping and loss rate are negligible for the 
short time used to determine the frequency. It is found experimentally that 
both the damping rate and the loss rate increase as Involve approaches the 
resonant value. 

In Fig. 26 (c) and (d) the results for the frequency as a function of involve 
are presented. The solid line shows the result of a two-body coupled-channels 
calculation of the molecular binding energy [80] . The parameters of the inter- 
atomic potentials are fit to the experimental results for the frequency. Clearly, 
the frequency of the coherent atom-molecule oscillations agrees very well with 
the molecular binding energy in vacuum over a large range of the magnetic 
field. Moreover, in the magnetic-field range involve — 157—159 G the frequency 
of the oscillations is well described by the formula |e m (.B)| = h 2 /ma 2 (B) for 
the binding energy, derived in Section 4.1.2. Close to resonance, however, the 
measured frequency deviates from the two-body result. The deviating exper- 
imental points are shown by open circles and are not taken into account in 
the determination of the interatomic potential. This deviation is due to many- 
body effects [81]. 

Although some of the physics of these coherent atom-molecule oscillations can 
roughly be understood by a simple two-level picture, it is worth noting that 
the physics of a Feshbach resonance is much richer. First of all, during Rabi 
oscillations in a simple two-level system one quantum in a state oscillates 
to the other state. In the case of a Feshbach resonance pairs of atoms oscil- 
late back and forth between the dressed-molecular condensate and the atomic 
condensate. Therefore, the hamiltonian is not quadratic in the annihilation 
and creation operators and the physics is more complicated. In particular the 
dressed molecule may decay into two noncondensed atoms instead of forming 
two condensate atoms. This process is discuss in detail below. Second, the 
observed atom-molecule oscillations are oscillations between an atomic con- 
densate and a dressed molecular condensate. The fact that one of the levels is 
a dressed molecule implies that by changing the magnetic field not only the 
detuning is altered, but also the internal state of the molecule itself. 

This is seen most easily by considering the linearized version of the time- 
dependent mean- field equation in Eq. (199). Writing a (t) = a e~^* /?i + 50 a (t) 
and 4> m (t) = <f) m e~ 2lflt / h + 5(f> m (t), we have that 



fl<tym(t) 
dt 




i\ ih- - 2/ffiHF fyjfi + 2g<PM,{t) , 
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(229) 



where we neglected the off-resonant part of the interatomic interactions. This 
is justified sufficiently close to resonance, where we are also allowed to neglect 
the energy-dependence of the atom- molecule coupling constant. 

Consider first the situation that the fractional derivative is absent in the lin- 
earized mean- field equations in Eq. (229), i.e., we are dealing with the model of 
Drummond et. al. [70], and Timmermans et al. [71,83]. These coupled equa- 
tions describe exactly the same Rabi oscillations as the coupled harmonic 
oscillators in Eq. (222), with the coupling equal to A = |4g0 a |. In the context 
of particle-number oscillations between condensates, Rabi oscillations are re- 
ferred to as Josephson oscillations and the associated frequency is called the 
Josephson frequency. The Josephson frequency in the absence of the fractional 
derivative term in Eq. (229) is given by 



hw b r = yj6 2 (B) + 16s 2 n a , (230) 

which reduces to fta>j arc ~ \5(B)\ sufficiently far off resonance where the cou- 
pling may be neglected. This result does not agree with the experimental 
result because, by neglecting the fractional derivative, which corresponds to 
the molecular self-energy, we are describing Josephson oscillations between 
an atomic condensate and a condensate of bare molecules instead of dressed 
molecules. Furthermore, using the result in Eq. (223) we have that the ampli- 
tude of these oscillations is given by 

/.bare _ . . 

Ai ~WW' { } 

In first approximation we take the dressing of the molecules into account as 
follows. If we are in the magnetic-field range where the Josephson frequency 
deviates not too much from the molecular binding energy, we are allowed 
to expand the propagator of the molecules around the pole at the bound- 
state energy. As we have seen in Section 4.1.3 this corresponds to introducing 
the dressed molecular field and leads to the Heisenberg equations of motion 
in Eq. (148). The linearized mean-field equations that describe the Josephson 
oscillations of a atomic and a dressed-molecular condensate are therefore given 
by 



. fc d6<j> m (t) 



ih ^ w =e m (£)<ty m (t) + 2g^Z{B)<t>M,{t) , 
»(*) 



it£?^=2gyfz(B)<j>*M m (t) , (232) 
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and lead to the Josephson frequency 



fiw 3 = + l6g 2 Z(B)n a , (233) 

which reduces to hujj ~ |e m (5)| in the situation where the coupling is much 
smaller than the binding energy. This result agrees with the experimental fact 
that the measured frequency is, sufficiently far from resonance, equal to the 
molecular binding energy. Moreover, the initial deviation from the two-body 
result in the measured frequency is approximately described by the equation 
for the Josephson frequency in Eq. (233). The amplitude of the oscillations is 
in this case given by 

16g 2 Z(B)n a . . 

A > - i^r • (234) 



which close to resonance is much larger than the result in Eq. (231). 

To get more quantitative understanding of the magnetic-field dependence of 
the Josephson frequency over the entire experimentally investigated range of 
magnetic field we calculate this frequency in a linear-response approximation, 
including the energy-dependence of the atom-molecule coupling and the atom- 
atom interactions. 

Before doing so, we make some remarks about the origin of the damping of 
the coherent atom-molecule oscillations and the overall loss of atoms that is 
observed in the experiments. One contribution to the damping is expected to 
be due to rogue dissociation [78]. Physically, this process corresponds to a pair 
of condensate atoms forming a dressed condensate molecule that then breaks 
up into two noncondensed atoms instead of oscillating back to the atomic 
condensate. This process is incorporated into our theory by the imaginary part 
of the molecular self-energy. As explained in Section 4.1.3 in the derivation 
of the Heisenberg equations of motion in Eq. (148), that involve the dressed 
molecules, we have neglected such a process. It is, however, incorporated in 
the full solution of the mean-field equation in Eq. (199). In the last section of 
the section we present the results of numerical solutions of these equations. 

The overall loss of atoms from the atomic condensate is also partially due 
to the rogue-dissociation process. The experimental fact that a significant 
thermal component is formed during the double-pulse sequence supports this 
idea. Apart from this process, it may also be that conventional loss processes, 
such as dipolar decay and three-body recombination play a role. Although such 
processes are expected to become more important near a Feshbach resonance, 
they are, however, not included in our simulations since there is no detailed 
knowledge about the precise magnetic-field dependence near the resonance. 
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Fig. 27. Molecular binding energy in vacuum. The solid line shows the result of 
a calculation with r^ g = 185ao- The dashed line shows \e(B)\ = h 2 /ma 2 . The 
experimental points are taken from [80]. The dotted line shows the detuning |<5(-B)|. 

In principle, however, these loss processes could be straightforwardly included 
in our calculations, by adding the appropriate imaginary terms to the mean- 
field equations. Another possible mechanism is the loss of atoms due to elastic 
collisions, the so-called quantum evaporation process [41]. This process is also 
not included in our present calculations. 

6.2 Josephson frequency 

With the mean-field theory derived in the previous sections we now calcu- 
late the magnetic-field and density dependence of the Josephson frequency 
of the coherent atom-molecule oscillations, in a linear approximation. The 
only parameter that has not been determined yet is the effective range of the 
interatomic interactions rb g . All other parameters are known for 85 Rb. 

The effective range is determined by calculating the molecular binding energy 
in vacuum and comparing the result with the experimental data. We have 
seen that far off resonance the Josephson frequency is essentially equal to the 
molecular binding energy. Since the effect of a nonzero effective range only 
plays a role for large energies, and thus is important far off resonance, this 
comparison uniquely determines the effective range. As explained in detail in 
Section 4.1.2, the molecular binding energy is determined by solving for E in 
the equation 

E - 6(B) - KE£\E) = 0. (235) 

For 85 Rb the background scattering length is negative and the effective range 
turns out to be positive. The retarded molecular self-energy is therefore given 
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(236) 
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In Fig. 27 the result of the numerical solution of Eq. (235) is shown for 
r bg = 185a - Also shown in this figure are the experimental data points. 
Clearly, far off resonance there is good agreement between our results and 
the experimental data points. Therefore, we use this value for the effective 
range from now on in all our calculations. The absolute value of the detun- 
ing is shown by the dotted line, and deviates significantly from the binding 
energy. The dashed line in Fig. 27 indicates the formula |e m | = h 2 /ma 2 . As 
we have derived in Section 4.1.2 this formula should accurately describe the 
magnetic-field dependence of the binding energy close to resonance. Clearly, 
the solid line that indicates the result that includes the nonzero effective range 
becomes closer to the dashed line as we approach resonance. However, there 
is a significant range of magnetic field where we need to include the effective 
range in our calculations. Closer to the resonance, the experimental points 
start to deviate from the two-atom binding energy. This deviation is taken 
into account by considering many-body effects. Note, therefore, that the ex- 
pected oscillation frequency h 2 /ma 2 never leads to a quantitative agreement 
with experiment. 

As mentioned previously, we calculate the many-body effects on the frequency 
of the coherent atom-molecule oscillations in linear approximation. Therefore, 
we first need to determine the equilibrium around which to linearize. In detail, 
the equilibrium values of the atomic and molecular condensate wave func- 
tions are determined by solving the time-independent mean-field equations 
in Eq. (185) together with the equation for the Hartree-Fock self-energy in 
Eq. (186) at a fixed chemical potential /i. To compare with the experimental 
results it is more convenient to solve these equations at a fixed condensate 
density. The chemical potential is then determined from these equations. 

In Fig. 28 we show the result of this calculation for an atomic condensate 
density of n a = 2 x 10 12 cm~ 3 . The solid line shows the Hartree-Fock self- 
energy ft£ HF and the dashed line the chemical potential as a function of the 
magnetic field, both in units of the energy ATia(B)h 2 n ai /m. Note that far off 
resonance, where the energy dependence of the interaction may be neglected, 
we have that /i = 4ira(B)h 2 n a /m and fi,S HF = 2/i. This is the expected result. 
The inset of Fig. 28 shows the fraction of bare molecules |</> m | 2 / n a- Note that 
this fraction is always very small. This justifies neglecting the atom-molecule 
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Fig. 28. Hartree-Fock self-energy (solid line) and chemical potential (dashed line) 
as a function of the magnetic field for an atomic condensate density of n a = 2 x 10 12 
cm~ 3 . Both quantities are shown in units of 4ira(B)ti 2 n a /m. Far off resonance, 
where the energy dependence of the interactions can be safely neglected we have 
that hT, UF = 8ira(B)h 2 n a /m and \i = 47ra(B)ti 2 n a ,/m, as expected. The inset shows 
the fraction of bare molecules as a function of the magnetic field. 
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Fig. 29. The dispersion relation for the collective modes of an atom- molecule system 
for a condensate density of n a = 2 x 10 12 cm~ 3 at a magnetic field of B = 157 G. The 
momentum is measured in units of the inverse coherence length = y/l6Tra(B)n a . 
The upper branch corresponds to the gapless dispersion for phonons. The solid line 
is the result of the full calculation, the dashed line shows the Bogoliubov disper- 
sion for the scattering length a{B). The lower branch corresponds to the coherent 
atom-molecule oscillations. The solid line is the result of the full calculation whereas 
the dashed line shows the result with the same zero- momentum part, but with the 
momentum dependence determined by /i 2 k 2 /4m. 



and molecule-molecule interactions since from this figure we see that the mean- 
field energies associated with these interactions are at least three orders of 
magnitude smaller. A posteriori this observation justifies neglecting the effect 
of the presence of the molecular condensate on the atoms in the approach of 
Kdhler et al. [84]. 
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Fig. 30. Josephson frequency of coherent atom- molecule oscillations for various 
values of the condendate density. The solid lines are the results of calculations for 
nonzero condensate density. The different lines correspond from top to bottom to 
the decreasing condensate densities n a = 5 x 10 12 cm -3 , n a = 2 x 10 12 cm" 3 , and 
n a = 10 12 cm -3 . The dashed line corresponds to the molecular binding energy in 
vacuum, i.e., n a = 0. The experimental data points, taken from Ref. [80], are also 
shown. 

Since the coherent atom-molecule oscillations are a collective mode where 
the amplitude of the atomic and molecular condensate wave functions oscil- 
late out-of-phase, we study the collective modes of the system. As explained 
in detail in the previous section, the frequencies of the collective modes are 
determined by Eq. (195). This equation is solved numerically and yields a 
dispersion relation with two branches. 

The result of this calculation is shown in Fig. 29 for an atomic condensate den- 
sity of n a = 2 x 10 12 cm" 3 and a magnetic field of B = 157 G. The momentum 
is indicated in units of the inverse coherence length £ _1 = yi67ra(£>)n a . The 
upper branch corresponds to the gapless phonon excitations. For small mo- 
menta this branch has a linear momentum dependence. The upper dashed 
line indicates the Bogoliubov dispersion in Eq. (191) evaluated at the scat- 
tering length a(B). For small momentum the solid and the dashed line are 
almost identical. For larger momenta the numerically exact result is smaller, 
due to the energy-dependence of the interactions that effectively reduce the 
scattering length. 

The lower branch corresponds to the coherent atom-molecule oscillations and 
is gapped. The solid line indicates the result of the full calculations. For small 
momenta it is well described by 

frw k ~ -ftwj + e k /2 , (237) 
where u>j is the Josephson frequency. The dispersion resulting from this last 
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equation is shown in the lower part Fig. 29 by the dashed line. This mo- 
mentum dependence is to be expected since sufficiently far from resonance 
the atom-molecule oscillations reduce to a two-body excitation. The fact that 
the dispersion is negative is due to the fact that we are linearizing around a 
metastable situation with more atoms than molecules. Although this is the 
experimentally relevant situation, the true equilibrium situation for negative 
detuning corresponds to almost all atoms in the molecular state [83]. 

In Fig. 30 we present the results for the Josephson frequency as a function of 
the magnetic field, for different values of the condensate density. The solid lines 
in this figure show, from top to bottom, the results for an decreasing nonzero 
condensate density. The respective condensate densities are given by n a = 
5 x 10 12 cm" 3 , n a = 2x 10 12 cm" 3 , and n a = 10 12 cm" 3 . The dashed line shows 
the molecular binding energy in vacuum. The Josephson frequency reduces 
to the molecular binding energy for all values of the condensate density, in 
agreement with previous remarks. Nevertheless, sufficiently close to resonance 
there is a deviation from the two-body result due to many-body effects. This 
deviation becomes larger with increasing condensate density. 

In order to confront our results with the experimental data we have to real- 
ize that the experiments are performed in a magnetic trap. Taking only the 
ground states </> a (x) and <j) m (x) into account for both the atomic and the molec- 
ular condensates, respectively, this implies effectively that the atom-molecule 
coupling g is reduced by an overlap integral. Hence we define the effective 
homogeneous condensate density by means of n & = N a [J rfx0 2 (x)0 m (x)] 2 = 
16A/2A^ a m 3//2 z/ rv /z^/(1257r 3 /i 3 / 2 ), where iV a denotes the number of condensed 
atoms and v r and v z the radial and axial trapping frequencies, respectively. 
For the experiments of Claussen et al. we have that N a ~ 8000 during the os- 
cillations close to resonance as seen from Fig. 26, which results in an effective 
density of n a ~ 2 x 10 12 cm" 3 . This agrees also with the effective homogeneous 
density quoted by Claussen et al. [80]. The solid curve in Fig. 30 clearly shows 
an excellent agreement with the experimentally observed frequency for this 
density. 

It is important to note that there are two hidden assumptions in the above 
comparison. First, we have used that the dressed molecules are trapped in the 
same external potential as the atoms. This is not obvious because the bare 
molecular state involved in the Feshbach resonance is high-field seeking and 
therefore not trapped. However, Eq. (133) shows that near resonance almost 
all the amplitude of the dressed molecule is in the low-field seeking open 
channel and its magnetic moment is therefore almost equal to twice the atomic 
magnetic moment. Second, we have determined the frequency of the coherent 
atom- molecule oscillations in equilibrium. In contrast, the observed oscillations 
in the number of condensate atoms is clearly a nonequilibrium phenomenon. 
This is, however, expected not to play an important role because the Ramsey- 
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Fig. 31. Josephson frequency of coherent atom-molecule oscillations as a function of 
the condensate density, for fixed magnetic field. We have subtracted the molecular 
binding energy. 



pulse sequence is performed on such a fast time scale that the response of 
the condensate wave function can be neglected. By variationally solving the 
Gross-Pitaevskii equation for the atomic condensate wave function, we have 
explicitly checked that after a typical pulse sequence its width is only a few 
percent larger than the harmonic oscillator ground state. 

Finally, we calculate the Josephson frequency as a function of the condensate 
density. The results of this calculation are presented in Fig. 31, for various 
values of the magnetic field which is kept fixed in these calculations. In the 
presentation of the results we have subtracted the molecular binding energy 
to bring out the many-body effects more clearly. As expected, the difference 
between the Josephson frequency and the molecular binding energy increases 
with increasing condensate density. Moreover, for values of the magnetic field 
closer to resonance the difference is also larger. 

The above calculations in the linear approximation give already a great deal of 
insight in the coherent atom-molecule oscillations, and, in particular, in their 
many-body aspects. In the next section we aim at achieving also insight in the 
nonlinear dynamics and damping resulting from the time-dependent mean- 
field equations for the double-pulse experiments. In particular, we also discuss 
the rogue-dissociation process. The nonlinear effects in these experiments has 
first been discussed by Kokkelmans and Holland [79], Mackie et al. [78], and 
Kohler et al. [84], on the basis of their mean-field approaches summarized in 
Section 5.2. 
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Fig. 32. Fraction of atoms in the atomic condensate. The solid line shows the result 
of the inclusion of the rogue-dissociation process into the calculations. The dashed 
line shows the result of a calculation without this process. The dotted line shows 
the result for a calculation that includes the estimate in Eq. (242). We have taken 
the parameters B m i t = 162 G, -B C voivc = 158 G, and n a = 2 x 10 12 cm~ 3 . 

6.3 Beyond linear response 

In this section we discuss the numerical solution of the time-dependent mean- 
field equations using the methods described in Sec. 5.1.2. We focus here on 
the situation where the detuning is only changed instantaneously, so that we 
are allowed to use the Green's function method discussed in this section. After 
the elimination of the molecular condensate wave function from the mean-field 
equations, the effective equation for the atomic condensate wave function is 
then given by 



In this equation, the term that involves the integral over frequencies describes 
the fact that a pair of condensate atoms that forms a molecule can decay 
into a pair of noncondensed atoms with opposite momenta, i.e., the rogue- 
dissociation process. In the absence of this term the equation effectively takes 
into account the dressing of molecules in an adiabatic manner, and describes 
Josephson oscillations between a condensate of atoms and dressed molecules. 



<90a(t) _ 47ra bg ^ 2 

dt m 



\Mt)\ 2 Mt) 




(238) 
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As we have discussed in the previous section, the above equation is only ap- 
plicable to the situation of a sudden change in magnetic field. Therefore, we 
perform the following calculation. For a given magnetic field B init and atomic 
condensate density we calculate the equilibrium values of the molecular wave 
functions and the Hartree-Fock self-energy, using the time-independent mean- 
field equations in Eq. (185) and Eq. (186). Then we change the magnetic field 
instantaneously to the value -Bevoivc and keep it at this value. In Fig. 32 the 
results of the calculations for this situation are shown, with B init = 162 G and 
-Bcvoivc = 158 G. The atomic condensate density is taken equal to n a = 2 x 10 12 
cm -3 . The dashed line shows the result for a calculation without the rogue- 
dissociation process and shows oscillations where a fraction of the atoms is 
converted into molecules and oscillates back and forth between the atomic 
and dressed molecular condensate. Since there is no decay mechanism, all of 
the atoms come back into the atomic condensate at times equal to a multiple 
of the oscillation period. The solid line shows the result of a calculation that in- 
cludes the rogue-dissociation process. Clearly, the number of condensate atoms 
oscillates in this well. However, not all of the atoms come back into the 

atomic condensate and there is a decay of the number of atoms in the atomic 
condensate. This is precisely due to the above-mentioned rogue-dissociation 
process. 

Although the preliminary calculations presented in this section are limited to 
the case of a step in the magnetic field, they nevertheless present some insight 
in the effects of the rogue- dissociation process on the coherent atom-molecule 
oscillations in a Ramsey experiment. In future work we intend to study also 
the case of time-dependent magnetic fields, by an exact numerical treatment 
of the fractional derivative in our time-dependent mean-field equations. In 
particular, we are interested in the magnetic-field dependence of the damping 
that is caused by the rogue-dissociation process. 

We can estimate this dependence as follows. The Green's function associated 
with the rogue-dissociation process, 



is sharply peaked in time. Hence we approximate this Green's function by 




2tt [hw + 2/l£ HF - 5(B)] 2 + (# 4 m 3 /47r 2 /i 6 )^ ' 



(239) 



o 



GM(t - t) ~ r(B)G<g(0)6(t - If) , 



(240) 
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with the timescale t(B) given by 

r{B)= f dt G«(t) , (241) 



with t c a positive cut-off that is determined such that the result for r(B) 
depends only very weakly on t c . The Green's function evaluated at zero time 
equals G^(0) = 1 — Z(B), a result which follows from the sum rule for the 
molecular density of states in Eq. (132). This gives the contribution 

_ -2i[l-Z(B)]g*T(B)^ ... 2 



\MWMt) , (242) 



to the right-hand side of Eq. (238). The rate equation for the atomic density 
that follows from this term is given by 

£ ^ _ 4,i - z(mMB) nl(t) (243) 

which after linearization leads to the following equation for the number of 
condensate atoms 

^ - -WAV • (244) 



with the rate j3 given by 

n 8[1- Z(B)]g 2 r(B)n a . . 

13 ~ V J± y V ; . 245 

nr 



We observe from this equation that the loss rate of atoms from the atomic 
condensate due to the rogue-dissociation process increases as the magnetic 
field approaches its resonant value. This is indeed what is observed experi- 
mentally [80]. Far off resonance the loss rate vanishes since the wave func- 
tion renormalization factor Z(B) — > 1 in this limit. For the parameters of 
Fig. 26 (a) at the effective homogeneous density n a = 2 x 10 12 cm -3 , we have 
that t{B) ~ 1.28 x 1(T 9 s, which leads to f3 ~ 0.45 kHz. The dotted line in 
Fig. 32 shows the result of a calculation that includes the term in Eq. (242). 
The exact result, shown by the solid line, and this approximate result show the 
same overall damping rate. This justifies the approximation for the Green's 
function in Eq. (240). The result for the damping rate (3 is about a factor of 
eight smaller than the experimental result. 
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Fig. 33. Fraction of atoms in the atomic condensate after a step in the magnetic 
field. The solid line corresponds to involve = 156.1 G. The dashed and dotted 
line correspond to a magnetic field of -B C voivc = 156.5 G and B cvo \ vc = 156.9 G, 
respectively. The initial magnetic field is -EWt = 162 G and the density of the 
atomic condensate is n a = 2 x 10 12 cm -3 . 

To further investigate the magnetic-field dependence of the damping of the 
coherent atom-molecule oscillations, we have calculated the numerical solution 
of the effective equation of motion for the atomic condensate wave function 
for a step in the magnetic field, for three different final magnetic fields. The 
results of these calculations are shown in Fig. 33. The solid, dashed, and dotted 
lines corresponds to a magnetic field of -B ev oive — 156.1 G, B cvo \ vc = 156.5 G, 
and -Bevoive = 156.9 G, respectively. The initial equilibrium corresponds to 
an atomic condensate density of n a = 2 x 10 12 cm" 3 at a magnetic field of 
-Binit = 162 G. Note the increase in the frequency with increasing magnetic 
field. 

The magnetic-field dependence of the frequency and damping of the coherent 
atom-molecule oscillations is found from these numerical results by fitting 
with the equation in Eq. (228). The results are presented in Fig. 34. The solid 
line corresponds to the Josephson frequency of the coherent atom-molecule 
oscillations that was found by means of the linear-response calculation of the 
previous section. The deviation for large magnetic fields is understood because 
we have, in our numerical solution of the effective mean-field equation, not 
taken into account the higher-order energy-dependences of the molecular self- 
energy that are fully taken into account in the linear-response theory. The 
inset shows the damping as a function of the magnetic field. Note the increase 
of the damping as the magnetic field approaches its resonant value. This is 
expected from the estimate in Eq. (245). 

The above analysis indicates that the rogue-dissociation process gives possibly 
a contribution to the experimentally observed damping of the coherent atom- 
molecule oscillations. Presumably, however, also other mechanisms contribute 
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Fig. 34. Frequency and damping as a function of the magnetic field. The solid line 
corresponds to the frequency found by means of linear-response theory. 

to the observed damping. In particular, we mention here the quantum evapo- 
ration process, that was shown to be important in the single-pulse experiments 
[41]. The detailed investigation of the damping of the coherent atom-molecule 
oscillation is a subject for further study. 



7 Conclusions and outlook 

In this review paper we have presented the derivation of an effective quan- 
tum field theory suitable for the description of a Bose gas near a Feshbach 
resonance, since it incorporates the two-atom physics exactly. We have pre- 
sented several applications of this theory, both above and below the critical 
temperature for Bose-Einstein condensation. In the last part of this paper we 
have studied in detail the magnetic-field dependence of the frequency of the 
coherent atom-molecule oscillations and have obtained excellent agreement 
with the experimental results. In particular, we have been able to quantita- 
tively explain the many-body effects on this frequency by making use of a 
linear-response approximation to our mean-field equations. Although we have 
already presented some numerical solutions of the mean-field equations that 
improve on this approximation, a great deal of work still has to be done. The 
numerical solution of these equations for the situation of time-dependent de- 
tuning is rather involved. Nevertheless, work in this direction is in progress 
and will be reported in a future publication. 

As already mentioned, we have also discussed the properties of the gas above 
the critical temperature. This discussion was mainly concerned with the equi- 
librium properties of the gas and we studied the many-body effects on the 
bound-state energy of the molecular state. An important conclusion of this 
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study is that, for certain values of the parameters, there exists a many-body 
induced resonant state with a relatively small energy. In future work we intend 
to study the effects of the appearance of this resonant state in the molecular 
density of states on the properties of the gas. In particular we expect that due 
to this effect the number of molecules in the gas will be large even at relatively 
large detuning, which can not be explained on the basis of two-atom physics. 

Furthermore, to study the normal state also in an out-of-equilibrium situation, 
we should derive a quantum kinetic theory that describes the evolution of the 
local occupation numbers of the atoms and molecules. Moreover, the descrip- 
tion of the Bose-Einstein condensed phase of the gas at nonzero temperatures 
requires a modification of the mean-field equations such that they include the 
effects of the thermal clouds of atoms and molecules, and we need equations 
for the evolution of the local occupation numbers of the latter. The extension 
of the theory presented in this paper to these situations can be derived in a 
unifying manner by using a functional formulation of the Schwinger-Keldysh 
nonequilibrium theory [113], and is especially important in view of the ongo- 
ing effort to produce ultracold molecules by means of a sweep in the magnetic 
field through the Feshbach resonance [62]. 

The theory presented in this paper is generalized to a gas of fermionic atoms 
in a straightforward manner [68,85]. One modification is that to have s-wave 
scattering between fermionic atoms we have to have a mixture of atoms with 
two hyperfine states, since the Pauli principle forbids s-wave scattering be- 
tween identical fermions. Furthermore, the properties of the dressed molecular 
state is altered due to the presence of the Fermi sphere. A molecule with zero 
momentum only decays if its energy is above twice the Fermi energy. If the 
molecular state lies below twice the Fermi energy, the equilibrium situation is 
a Bose-Einstein condensate of molecules. If we start from this situation and 
increase the detuning, the Bose-Einstein condensate of molecules crosses over 
to a Bose-Einstein condensate of Cooper pairs, i.e., a BCS-BEC crossover oc- 
curs [61,57]. In view of the ongoing experiments with atomic Fermi gases near 
a Feshbach resonance [62,63,64,65,66,67], it is particularly interesting to study 
the effects of nonadiabticity on the crossover from a Bose-Einstein condensate 
of molecules to a degenerate Fermi gas. In particular, the atomic distribution 
function after such a sweep, and its dependence on the duration of the sweep, 
is of great interest, since this will determine whether or not a BCS-state will 
form after equilibration. Determination of the atomic distribution function 
requires, in first instance, knowledge of the solution of the mean-field equa- 
tion for the molecular condensate for time- dependent detuning. Work in this 
direction is in progress. We also intend to study the equilibrium properties of 
this crossover, and in particular the behaviour of the critical temperature, in 
detail in future work. 

Clearly, Feshbach resonances present an exciting opportunity for the exper- 
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imental and theoretical study of the many-body properties of atomic and 
molecular Bose and Fermi gases. There is little doubt that these Feshbach 
resonances will find many new applications in the years to come. 
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